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ABSTRACT 

We present the Strong Lensing Legacy Survey - ARCS (SARCS) sample compiled from the final 
T0006 data release of the Canada- France-Hawaii Telescope Legacy Survey (CFHTLS) covering a total 
non-overlapping area of 159 deg^. We adopt a semi-automatic method to find gravitational arcs in the 
survey that makes use of an arc-finding algorithm. The candidate list is pruned by visual inspection 
and ranking to form the final SARCS sample. This list also includes some serendipitously discovered 
lens candidates which the automated algorithm did not detect. The SARCS sample consists of 127 
lens candidates which span arc radii ~ 2" — IS" within the unmasked area of ~150 deg^. Within 
the sample, 54 systems are promising lenses amongst which, we find 12 giant arcs (length-to- width 
ratio > 8). We also find 2 radial arc candidates in SL2SJ141447+544704. From our sample, we detect 
a systematic alignment of the giant arcs with the major axis of the baryonic component of the putative 
lens in concordance with previous studies. This alignment is also observed for all arcs in the sample 
and does not vary significantly with increasing arc radius. The mean values of the photometric 
redshift distributions of lenses corresponding to the giant arcs and all arcs sample are at z ^ 0.6. 
Owing to the large area and depth of the CFHTLS, we find the largest sample of lenses probing mass 
scales that are intermediate to cluster and galaxy lenses for the first time. We compare the observed 
image separation distribution (ISD) of our arcs with theoretical models. A two-component density 
profile for the lenses which accounts for both the central galaxy and the dark matter component is 
required by the data to explain the observed ISD. Unfortunately, current levels of uncertainties and 
degeneracies accommodate models both with and without adiabatic contraction. We also show the 
effects of changing parameters of the model that predict the ISD and that a larger lens sample might 
constrain relations such as the concentration-mass relation, mass-luminosity relation and the faint-end 
slope of the luminosity function. 

Subject headings: dark matter - gravitational lensing: strong - methods: data analysis - surveys: 
CFHTLS-SL2S 



1. INTRODUCTION 

Gravitational lensing is the deflection of light com- 
ing from distant sources in the Universe, due to the 
gravitatio nal potential of intervening s t ructu r es (see re- 
views e.g. IBlandford NaravanI [1981 [1991 IKochanekl 
|2006[ ). The last decade has seen the rise of a wide 
variety of applications of strong lensing such as the 
study of distant lensed galaxies with unprecedented mag- 
nification (e.g., Jmpellizzeri et al.i 12008; .Swinbank et all 
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[2QQ1 iZitrin fc BroadhurstI [2QQ1 iRichard et all [Mil ), 
the constr aints on s u bstruct u re within lensing ha- 
los (e.g., 'More et al.' I2009al : ISuvu fc Halkolal 120101 : 
IVegetti et al. 2010a. b), accurate ni e asurements of th e 
Hubble constant (e.g., IColesI 120081 : ISuvu et al.l I2010D . 
constraints on the stellar initial mass function (e.g., 
Treu et al.l [20101 : iFerreras et all 120101 : ISonnenfeld et al.l 
201 1[ ). constraints on the slope of the inner density 
profile of the lensing halos (e . g., iTreu fc KoopmansI 



2002a, b: Koopmans fc Treu 2003; Koopmans et al. 2006 
More e t al. 2008; Barnabe e t al. 2009; Koopmans et al.l 
2009h and estimation of the fraction of dark mat- 
ter in gal axy-scale halos (e.g., 'Gav azzi et al.l 120071: 
Jiang fc Ko chanek 2007; Grillo et al. 20101: iFaure et al.l 



20111 : IMore et al. h2011ai : iRuff et al.li201l V 



Although strong lensing is a rare event, several sur- 
veys covering a wide sky area and deep enough imag- 
ing across different wavelengths, have resulted in the 
discovery of over 200 strong lens systems at galaxy 
scales from surveys such as, the Great Observatories Ori- 
gins Deep Survey (Fassnacht et al. 2004), Cosmic Evolu- 
tion Survey (COSMOS Faure et al. 2008; Jacksonil20^, 
Mediu Deep Survey (Ratnatunga et al. 199 9) the Cos- 
mic Lens Ah Sky Survey (CLASS, iMvers e t al. 2 QQ3I) 
and the Sloan Lens ACS Survey (SLAGS, Bolt on et al.l 
|2006[ ) and about a few dozen lens systems at cluster 
scales such as, the Extended Medium Sensitivity Sur- 
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vey fEMSS. [LuDDino et al J'lQQQ"). the MAssive Cluster 
Survey (MACS. lEbeling et al. 2001 ), the Las Campanas 
Distant Cluster Survey (LCDCS, IZaritsky fc Gonzalezl 
12003) an d the Red sequence Cluster Survey (RCS, 
[Cladders et alj l2003l . henceforth G03). Large imag- 
ing and spectroscopic surveys enable us to probe sta- 
tistical properties of both dark matter and bary- 
onic matter or constrain cosmological parameters with 
high accuracy. For instance, on galaxy scales, the 
SLACS sample has been used to study the aver- 
age density profiles of lens galaxies up to redshift of 
■^0.3 (e.g., [Koopmans et al.l 120061: IGavazzi et al.l 120071: 



iKoopmans et al.l 120091 : lAuger et al.l l2010f ). iFalco et all 
(| 19991 ) used the Cf A- Arizona Space Telescope LEns Sur- 
vey (CASTLES, iMuhoz et all 11998 ^ sample to measure 
the mean exti nction due to the interstellar medium of the 
lens galaxies, Mediavilla et al. (2009) estimated the frac- 
tion of mass in compact objects within lens galaxies from 
the CASTLES, IWxithd (|20"04D constrained the bright end 
of the quasar luminosity function from the absence of 
lensed quasars at high redshift in the Sloan Digital Sky 
Survey (SDSS). On cluster scales, G03 found the ob- 
served abundance of giant arcs from the RCS is too high 
to be consistent with the predictions from th e current 
standard cosniolo^ical mo del f see also .Bartelmann et al.l 
IT998I : iLi et aLll2005l l2006l) and lZitrin etHI (l20TTI) found 
some discrepancy between the predictions from the stan- 
dard model and the observed distribution of Einstein 
radii from the MACS sample. However, the magnitude 
of these diff erences has been ni i tigated with subsequen t 
studies fe.g.. lHoresh et al.ll2005l : iMeneghetti et al.ll201ll ). 

As discussed above, the majority of the surveys in the 
past have primarily focused on studying galaxy-scale or 
cluster-scale structures. As a result, matter distribution 
in galaxies and galaxy clusters is relatively well-studied 
via both strong and weak lensing. A further improve- 
ment in our understanding has come from the use of com- 
plementary methods to lensing such as stellar kinematics, 
satellite kinematics and X-ray scaling relations. In con- 
trast, little is known about galaxy groups which are in- 
termediate to galaxies and galaxy clusters, typically cor- 
responding to masses of 10^2 - 10^"^ Mq. Relatively fewer 
investigations have been carried out with galaxy groups 
e.g., study of intra-^ roup medium with v ery low redshift 
X-ray sample (Helsd on fc Ponmanll2000[ ). study of mass- 
to-light ratios with the Ca nadian Network for Observa- 
tional Cosmology 2 sample (|Parker et al.l l20Q5[). study of 
faint end of the luminosity function of nearby compact 
groups (;Krusch et al.ll2006f ). st udy of concentration-mass 
(c-M) relation from the SDSS (|Mandelbaum et al.ll2QQ8l ) 
via weak lensing, study of colors and star formation (e.g., 
[Balogh et al. 2009, 2011), study of scaling relations of X- 
ray selected groups (Rines fc Diaferio 201 Qj) and study 
of baryon fractions from the 2MASS (D ai et al.l 12010). 
Since studies of groups are limited and we still do not 
have a detailed understanding of matter distribution, for- 
mation and evolution of galaxy groups. Being one of 
the important components in the hierarchical assembly 
of structures in the Universe, galaxy groups are much 
more massive than galaxy-scale halos and are concen- 
trated enough to act as lenses. Furthermore, since galaxy 
groups are quite abundant compared to massive struc- 
tures like galaxy clusters, the probability to find group 
scale lenses is also large. Hence, lensing can be success- 



fully used to study group-scale halos. 

The Stron g Lensing Legacy Survey (SL2S, 

iCabanac et al.l 1200 7) is a survey from the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS). 
The design of the CFHTLS allows us to find large 
sample of group-scale lenses, which can be studied in 
detail upto high redshifts, for the first time. The SL2S is 
as a precursor to wide field imaging surveys such as the 
Dark Energy Survey, Large Synoptic Survey Telescope 
and Euclid. A combined study of SL2S galaxy-scale 
lenses with the SLACS sample have been used to show 
mild evolut ion of the slope o f the average density profile 
of galaxies (iRuff et al.ll201lt ) which is constrained iDy the 
strong lensing and stellar dynamics techniques. Disk 
galaxies are also a relatively less studied population 
especially at high redshifts. An automated search 
for edge-on disk lenses from the SL2S resulted in 18 
candidates , out of which 3-5 are expected to be real 
lenses (Svg net et al]|2010f ). A subset of the SL2S groups 
have been studied in detail using a combination of 
techniques such as strong lensing, group dynamics, 
weak lensing to probe t he dens i ty pro fi les of the lensing 
groups (e.g., L imou sin et al.l 120091 : iThanjavur et al l 
12011 iVerdugoet al.ll2QTT[ ). 

In this paper, we present the SL2S-ARCS sample from 
the final T0006 release of the CFHTLS. In Section 2, we 
give an overview of the survey and procedure of sample 
selectio n, describe details of the algorithm, arcfinder 
("Alard" 2006) and present the final sample. In Section 3, 
we discuss some statistical results using the final sample. 
In Section 4, we summarize the survey and our main 
findings. 

2. THE CFHTLS-SL2S ARCS SAMPLE 

In this section, we give a brief overview of the survey 
from which we derive the lens sample. This is followed by 
a description of the semi-automatic process of selecting 
the candidates in the final sample. We also discuss how 
the algorithm, ARCFINDER works and the modifications 
implemented in the new arcfinder. Lastly, we present 
the final sample and report duplicate detections of some 
candidates from other surveys. 

2.1. Survey Overview 

CFHTLS is a photometric survey made with the 
Canada-France-Hawaii Telescope (CFHT) in five optical 
bands {u* g' r' i' z') using the wide-field imager MegaPrime 
with a field-of-view of 1 deg^ on the sky and a pixel 
size of 0.186''. The WIDE and DEEP components 
of the CFHTLS are designed to carry out extragalac- 
tic research. These components are ideal for searching 
strong lens systems. The SL2S sample is compiled from 
the CFHTLS-WIDE encompassing a combined area of 
171 deg^ and CFHTLS-DEEP encompassing a com- 
bined area of 4 deg^. However, taking into account the 
masked and overlapping areas, the effective area of the 
survey is 150.4 deg^ (146.9 deg^ for WIDE and 3.5 deg^ 
for DEEP). The WIDE consists of four fields Wl, W2, 
W3 and W4. The field Wl has the largest sky coverage of 
63.65 deg^. The fields W2 and W4 have similar sky cov- 
erages of 20.32 deg^ and 20.02 deg^ respectivel^fl The 



^ These numbers are est imated 

|http: / / terapix.iap.fr / cplt / table_syn_T0006.html| 
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field W3 has a sky coverage of 42.87 deg^ and is more 
than twice as large as W2 and W4. The DEEP also 
consists of four fields Dl (located within the Wl field), 
D2, D3 and D4. Each of the deep fields covers an area 
of 1 deg^. The DEEP images are produced in two image 
stacks D-25 and D-85. The former consists of 25% of 
the best seeing images and the latter consists of 85% of 
best seeing images. We use the D-25 images for search- 
ing lens candidates. Among the WIDE fields, g band 
imaging is the deepest of all bands with a limiting mag- 
nitude of 25.47 and a mean seeing of 0.78'' whereas the g 
band imaging of DEEP fields has nearly 10 times deeper 
exposures than the WIDE fields and the median seeing 
is ~ 0.7''. The zero point to convert fiux to AB mag- 
nitude for all bands is 30. Further details of the T0006 
release, which is the first complete release of the WIDE 
and DEEP, can be found on Terapix websit43. 

2.2. Sample selection 

The SL2S lens sample is compiled using two algo- 
rithms: RINGFINDER and ARCFINDER. The former aims 
at detecting galaxy-scale lenses by using color informa- 
tion. The SL2S RINGS sample will be presented and 
discussed in a separate paper (Gavazzi et al., in prepa- 
ration). Here, we focus on the SL2S-ARCS (SARCS) 
sample created with the help of arcfinder, followed by 
visual inspection and screening of the candidates. 

We define the SARCS sample such that lens systems 
with arc radius {Ra) ^ 2" belong to the sample. The 
radius of the arc is defined as the distance of the lensed 
image from the putative lens galaxy which is roughly 
the Einstein radius. Typically, lensing halos with Ein- 
stein radius larger than 2"are very massive lenses with 
significant contribution from the environment of the pri- 
mary lensing galaxy. Thus, the SARCS sample, pre- 
dominantly, consists of group to cluster scale lenses. Lens 
systems with Ra < 2", typically, form part of the RINGS 
sample. We note that a few lens systems are common to 
both SL2S samples since the cut on Ra is not a sharp 
limit imposed by the algorithms. 

The SARCS sample from the CFHTLS is compiled in 
a three-step process. The first step is to run the arc- 
detection algorithm called arcfinder. We choose to 
run the ARCFINDER on the ^-band image since most of 
the lensed images correspond to galaxies with high star 
formation that have little emission at redder wavelengths. 
Also, focusing on g-band prevents from detecting high 
redshift g-dropouts. However, we plan to search for g- 
dropout arcs which are brighter in the i band. These 
results will be presented in a separate paper. 

At the end of the first step, we produce a list of arc 
candidates with various parameters. The next step is ap- 
plying a cut-off on arc properties such as the area, the 
peak fiux count, the radius of curvature and reject candi- 
dates within masked area. These cuts allow a significant 
reduction in false detections at the cost of losing some 
real arc candidates. In the third step, visual inspection 
and classification are carried out to grade the quality of 
the candidates. Note that the final SARCS sample con- 
sists of candidates that are detected by the arcfinder 
and/or by visual inspection. 

^ |http://terapix.iap.fr/cplt/T0006-doc.pdf| 
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Fig. 1. — Flowchart showing the arcfinder algorithm. 



2.3. Automatic detection of arcs: arcfinder 

With the advent of several large imaging surveys, 
searching for lens systems visually is a subjective and 
time expensive exercise which may not be easily repeat- 
able. Hence, several algorithms have been devised in the 
last few years to automate the process of lens detection 
as much as possible. The algorithms, usually, focus on 
a specific target population, for example, algorithms to 
detect the lensed sourc es such a s quasars via the quasar 
time variability (Koch anek et aL 2006) based on the dif- 
ference imaging technique of lAlard fc LuptonI ()l998f ). 
spectrum-bas ed algorithms such as the one used by 
iBolton et all (|2004) to form the SLAG S sample following 
the technique of Warre rTet al.l (11996 *) or lens-modelling 
robots that assesses th e probabilit y of a bright red galaxy 
being a lens (Mars hall et al.ll200 9). Since groups to clus- 
ter scale lenses form arc-like lensed images, there are a 
few arc-finders in the literature fe.p!;.. iLenzen et al.|[2004l : 
iHoresh et aH 120051 : ISeidel fc Bartelmannll2007[ ) that aim 
at detecting elongat ed arc-like im ages. 

The ARCFINDER (Alard 2006) is a generic algorithm 
that can be used to detect elongated and curved fea- 
tures in an image. The algorithm uses pixel intensities 
from a standard FITS image to trace the structure of 
a feature. Multiple thresholds are applied to the struc- 
tural properties of the feature to select an arc candidate. 
The reader is referred to Alard (2006) for the details of 
the algorithm. Below, we describe the algorithm along 
with some modifications implemented in the newer ver- 
sion (V2.0) which is used to compile the SARCS sample. 

We show a fiowchart to illustrate the various steps in- 
volved in the algorithm (see Fig. [1]) and step through an 
example of a mock arc image (see Fig. [2]). 

• We run the arcfinder on images of 19354x19354 
pixels which corresponds to a single MEGACAM 
pointing. For such large images, the assumption of 
Gaussian distribution for the noise holds well. As- 
suming a Gaussian distribution for the noise in the 
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Fig. 2. — Outputs of arcfinder at various stages of the execution (panels 1-5) using a mock arc image. Panel 1 - original image with 
the mock arc, panel 2 - smoothed image, panel 3 - (disconnected) pixels with elongation estimator values higher than a threshold, panel 4 
- connected pixels above thresholds that belong to the arc and panel 5 - the reconstructed shape of the final arc. The image is 11.3^^on the 
side. Schematic diagram for measuring length of an arc (panel 6). 



input image (panel 1 of Fig. [2]), the noise a is cal- 
culated in the first step. Most of the arcs are unre- 
solved with the ground-based telescopes like CFHT 
and are limited by the size of the point spread func- 
tion (PSF). 

• Convolving the image with a smoothing kernel, 
with a size of the order of the PSF, can help damp 
the noise and enhance the detection of arcs. Hence, 
the input image is smoothed with a Mexican-hat 
filter in the next step (panel 2 of Fig. [2]). 

The Mexican-hat band-pass filter with scale b is 
given by 



M{x,y) = e ^ 
where = 3 pixels and 



0.5e- 



(1) 



(2) 



• Now, we define a local estimator of elongation for 
every pixel in the smoothed image. This is cal- 
culated by using the flux from all pixels within a 
window centered on every pixel. The size of the 
window is chosen to be a few times the PSF size. 
We consider a square window of side W centered on 
the pixel at (x,y). Using the second order moments 
of brightness distribution within the window, the 
local direction of the elongation of the feature is 
calculated. This direction is then used to align the 
feature along the x-axis to determine the local axis 
ratio, similar to a length-to-width ratio, within the 
window. The elongation estimator is defined as fol- 
lows 



TABLE 1 

Thresholds used in the final selection of the SARCS 
candidates. 



Step 


Parameters and Thresholds 


Smooth 


W=9, b'^ = 3.0 


Detect 


W=ll, Thl=ll*0.75(7 


Connect 


W=9, Th2=1.25 




W=9, Th4=0.5Sl,0.3S2 


Properties 


Final Thresholds 


Area 


55 < A < 500 (pix) 


peak counts 


< 50 (ADU s-i) 


mean counts 


2o- < SB < 50 (ADU s-^ pix-^) 


length 


> 7 (pix) 


width 


1.5 < w < 8 (pix) 


curvature 


rc < 100 and rc > 1000 (pix) 



W Fy^ 



(3) 



Note. — W is the window size in pixels. Thl, Th2 and Th4 
imply thresholds 1, 2 and 4, respectively. SI and S2 imply interme- 
diate mean counts of candidates in the input and smoothed images, 
respectively. 



where Fx is the integrated flux of the single central 
row of pixels and i^Y^ax maximum of the in- 

tegrated flux values calculated from single columns 
of pixels within the window. If Fx and Fy^^^ are 
above {W + l)0.75cr (threshold 1 in Fig. □ and Ta- 
ble [1]), then the estimator Q{x^ y) is assigned to the 
pixel (x,y) at the center of the window else the pixel 
is assigned a value of 0. Likewise, the estimator is 
calculated for every pixel in the image (see panel 3 
of Fig.E]). 

• In the following step, we attempt to connect the 
pixels that possibly belong to the candidate arc. 
A pixel is accepted if the estimator value from the 
image in panel 3 is > 1.25 (threshold 2 in Fig. [T] 
and Table [T]). The estimator is set to as soon as a 
pixel is connected in the prior step to avoid repet- 
itive checking of connected pixels. Iterating the 
above steps over the image allows tracing of the 
primary shape of the arc. If the number of con- 
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nected pixels is more than 10 (threshold 3) then 
the properties of arcs such as the length, peak flux 
and surface brightness are calculated. The panel 4 
shows the image with connected pixels which, for 
the case at hand, consists of one candidate only. 
The estimator is particularly suited for recovering 
pixels along the length of the arc. In order to con- 
nect the pixels along the width of the arc, we use 
the input and smoothed images (panels 1 and 2, re- 
spectively). Using arcs with length > 7 pixels, we 
apply surface brightness thresholds on pixel values 
from images in panels 1 and 2 in order to accept 
pixels belonging to the arc (threshold 4, see Fig. [T] 
and Table [T]) and construct the arc fully. 

• In the last step, if the candidate arc satisfies thresh- 
olds on the arc properties such as the width (> 1.5 
and <8 pixels), area (> 25 pixels), surface bright- 
ness (> 2(j) and peak fiux (< 500 ADU/s/pixel) 
which correspond to threshold 5 in Fig. [H then the 
arc is accepted (see panel 5 in Fig. [2j). 

• In the following, we describe how the arc proper- 
ties are measured. The area of the arc is the total 
number of pixels belonging to the final arc. The 
length of the arc is calculated by assuming that the 
candidate arc is similar to an arc of a circle. We 
first find the extreme ends of the arc, connect them 
with a chord labeled AB in the panel 6 of Fig. [51 
From the midpoint E of chord AB, we draw a line 
perpendicular to the chord and find its intersection 
D with the arc. The lengths ED and AB uniquely 
identify the circle to which the arc belongs and as 
well as the length of the arc. The width of the arc 
is defined as the ratio of the area to the length of 
the arc. The radius of the circle going through the 
arc is used as a proxy for the curvature of the arc. 

In the algorithm, the use of a second image at another 
wavelength and a mask file are optional. The algorithm, 
in its current version, merely prints the values of the 
pixels belonging to the candidate arc from the image 2 
and/or the mask file, if provided. This allows us to op- 
tionally screen the candidate arcs based on their color 
and/or masking information. For the SARCS sample, 
we made use of the mask option only and further re- 
stricted the parameter space of the output list of arcs by 
introducing more strict cuts on some of the arc proper- 
ties. The Table [T] lists the final set of thresholds which 
are satisfied by the SARCS candidates detected by AR- 

CFINDER. 

The components enclosed by the dashed line in the 
Fig. [T] are the sections of the algorithm that have been 
improved. This includes the way in which the pixels be- 
longing to the arc are connected and the way the proper- 
ties of the arcs are measured. In the earlier version of the 
algorithm, a candidate arc is accepted at the final stage 
only if it satisfies a certain threshold on the curvature of 
the arc. This feature has been removed and the option of 
using mask information is introduced in the V2.0. The 
values for all the different thresholds are tuned from an 
initial sample of CFHTLS arcs, found visually or from 
the previous version of arcfinder, based on an early 
release of the data. Furthermore, all the thresholds can 
now be set during the execution of the algorithm. 



The ARCFINDER V2.0 finds ~ 1.5 x more detections, 
most of which arise due to spikes and halos near stars. 
However, the use of masks removes these false detections 
thereby making the number of detections comparable or 
less than the output of the earlier version. The modified 
ARCFINDER is over 3 times faster than the earlier version. 
We have also recently parallelized the algorithm which 
enables us to achieve even faster computation on shared 
memory platforms. The ARCFINDER V2.0 is available 
upon request to the author. 

Finally, we note that the existing algorithms are far 
from perfect and almost always require manual interven- 
tion. Attempts to increase the completeness of the sam- 
ple almost always leads to a corresponding increase in 
the rate of false positives. In light of these issues, citizen 
science projects such as the Galaxy Zoo (jLintott et al.l 
2008) might be a better tool, for the time being, in identi- 
fying lens systems not detected by the algorithms. These 
projects, in turn, could be used to calibrate and improve 
existing algorithms. We are currently pursuing the fea- 
sibility of such a project. 

2.4. SARCS sample 

After the automatic detection and screening, about 
1000 candidates/deg^ are visually inspected. This is re- 
duced to a total sample of 413 candidates which is further 
considered for ranking by three people. The individually 
assigned ranks are from 1 to 4 with 1 being the least likely 
to 4 being the almost certain lens system. We present 
the 127 lens candidates which are ranked 2 or higher 
(average of the ranks by three people) and which have 
Ra > 2''. Ah of the 127 candidates E3 are listed in Ta- 
ble [21 that gives the ID, lens position, lens magnitudes in 
AB, lens redshifts with la uncertainties, arc radii, ranks, 
the type indicating whether the candidates are primarily 
detected via the arcfinder and the field name in which 
the candidate is located. For the calculation of arc radii 
in physical units, we use fiat cosmology with following 
parameters (1^^, I^a, i^o) = (0.3,0. 7,100 km s"^ Mpc"^). 
The symbol A stands for detection by arcfinder and V 
implies candidate is found serendipitously. A total of 54 
candidates are good to best systems (that is, rank of 3 
and above) and are shown in the Fig. [3] and Fig. HI Out 
of the total sample, 27 systems have been or are being 
followed up for further analyses (see Table [3]) and most 
of them are confirmed lens systems, 5 systems are very 
likely lenses and rest of the 96 are possible lenses. The 
Table [3] does not report information about any archival 
data on any of the lens candidates. 

One of t he deep fields , D2, has an overlap with the 
COSMOS fSco ville et al.l l2007). a photometric survey of 
1.8 deg^ in the / band (limiting magnitude of 26.5) with 
the Advanced Camera Survey (ACS) on HST. We find 
two of the SARC S candidates com mon to the lens sam- 
ple of COSMOS (jFaure et al.ll2008f ). The SARCS can- 
didates with IDs SA78 and SA83 are known COSMOS 
candidates and the rest of the 6 candidates within the 
COSMOS field are new (see Table [2] with Field labelled 
as D2). For comparison, we also show the ACS images 
of the SARCS candidates in Fig. [5] from the COSMOS 
archive. The W3 field overlaps with one of the 22 indi- 

High resolution images of these systems are made available 
at http:/ /kicp. uchicago.edu/~anupreeta/sarcs_sample 
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vidual fields of the RCS survey which covers a combined also detected in the SARCS sample identified by SA102 
area of ~90 deg^. The overlapping RCS field has a lens- (see Table [2]). 
ing cluster, RCS 1419.2+5326 at z = 0.64. This cluster is 



TABLE 2 Characteristics of lens candidates from the SARCS sample. 



ID 


RA 


Dec 


9 


r 


I 


^phot 






Rank 


Type 


Field 




hms 


dms 


mag 


mag 


mag 




/i-^Kpc 








C A 1 


no.m .01 on 
Uz:Ui:zi.oy 


nn. 1 . 1 cc nn 

— uy:io:io.uy 


22.06 


20.58 


19.90 


n A /^_Ln no 
U.4b±U.Uz 


2.2 


9.0 


2.0 


A 

A 


W 1 




no.no. 1 n ccn 
Uz:Uz: iU.oU 


1 1 .nn. 11 ci.Q 
— ll:Uy:ll.bo 


19.79 


18.48 


17.80 


n A o_Ln no 
U.4o±U.Uz 


5.0 


20.9 


Q 

0.7 


A 

A 


W 1 


C A Q 

oAd 


no.no. QQ QT 

Uz:Uz:c)o.o7 


n/^ .0 A .t^a 10 
— Ub:o4:ob.lz 


20.89 


19.91 


19.54 


n Q'T-Ln nQ 
U.o7±U.Uo 


2.2 


n 

7.y 


2.0 


A 

A 


W 1 


Q A A 

bA4 


no.nQ.no qa 
Uz:Uo:Uz.o4 


no. 01 .1/1 OK 
— Uo:zi: i4.zo 


01 C\Ci 

zi.yb 


on cc'T 
zU.O / 


1 n nn 

iy.yy 


n 1 /I J_n n'y 
U. i4zLU.U / 


A 

Z.4 


A 1 

4. i 


n 
z.U 


V 


W i 


c A c; 
oAo 


no. no. 10/^1 
Uz:Uo: iz.bi 


1 n. /I T.c^T nc^ 
— lU:47:U7.yo 


22.02 


20.55 


19.52 


n <^o_Ln nQ 
U.bz±U.Uo 


3.0 


1 /I Q 

14. 


2.3 


^ T 
V 


W 1 


Q A 
oAD 


UZ.Uo.ZU.4o 


07. Q/l . cn 70 
— U < .o4.0U. I 


01 <^Q 

zi.by 


on Q<^ 
zu.ob 


1 Q /I c; 

iy.4o 


n c;Q_|_n HQ 

u.oyztu.uo 


c; n 
O.U 


OQ 

zo.z 


Q Q 
0.0 


A 


VV i 




n0.nQ.4Q 


nQ.4o.c;Q c;i 


1 7 R4 


lU. ( 


1 1 R 
lu. 10 


n oc:_|_n no 
u.zozcu.uz 


^ 

O.U 


1 7 
10. ( 


0.0 




W1 
VV 1 


SA8 


02:04:54.51 


— 10:24:02.48 


19.52 


18.18 


17.68 


33+0 02 


10.8 


35.9 


2.7 


A 


Wl 


SA9 


02:05:03.15 


-11:05:46.63 


21.54 


20.17 


19.16 


0.62±0.03 


3.3 


15.7 


3.0 


V 


Wl 


SAIO 


02:06:48.47 


-06:57:01.33 


20.95 


19.60 


18.90 


0.49±0.06 


3.2 


13.5 


3.0 


A 


Wl 


SAll 


02:08:15.66 


-07:24:57.81 


21.98 


20.55 


19.49 


0.62±0.02 


4.3 


20.4 


2.0 


A 


Wl 


SA12 


02:08:16.87 


-09:36:52.69 


22.16 


20.85 


19.60 


0.74±0.03 


3.4 


17.4 


3.7 


A 


Wl 


SA13 


02:08:41.61 


-07:01:28.07 


19.85 


18.70 


18.20 


0.29±0.03 


3.5 


10.7 


2.0 


V 


Wl 


SA14 


02:09:29.33 


-06:43:11.26 


20.47 


19.15 


18.57 


0.45±0.02 


3.2 


12.9 


3.7 


A 


Wl 


SA15 


02:09:57.67 


-03:54:57.08 


21.41 


19.98 


19.27 


0.43±0.03 


3.9 


15.3 


2.3 


V 


Wl 


SA16 


02:10:26.57 


-04:46:41.59 


21.77 


20.79 


20.24 


0.55±0.03 


1.9 


8.5 


2.0 


A 


Wl 


SA17 


02:10:51.59 


-03:52:52.64 


21.79 


20.84 


19.91 


0.73±0.04 


1.9 


9.7 


2.0 


V 


Wl 


SA18 


02:11:08.66 


-10:12:13.86 


19.81 


18.44 


17.89 


0.38±0.02 


2.0 


7.3 


2.3 


A 


Wl 


SA19 


02:11:18.49 


-04:27:29.20 


23.13 


22.48 


21.43 


1.19±0.07 


3.5 


20.3 


3.3 


A 


Wl 


SA20 


02:12:20.52 


-09:38:44.10 


23.48 


22.18 


20.88 


0.77±0.03 


2.4 


12.4 


2.0 


A 


Wl 


SA21'^ 


02:13:24.52 


-07:43:54.82 


24.11 


23.73 


23.18 


0.80±0.16 


2.8 


14.7 


4.0 


V 


Wl 


SA22 


02:14:08.07 


-05:35:32.39 


20.98 


19.42 


18.69 


0.48±0.02 


7.1 


29.7 


4.0 


A 


Wl 


SA23^ 


09-1 4-1 1 94 

UZ . J.^. J. J. . z^ 


— n4-nt=;-n9 71 

— U^.UO.UZ. ( J. 


22.11 


20.92 


19.88 


n 7/1 -un c\A 
U. ( ^zcu. u^ 


1.9 


9.7 


4.0 




Wl 


SA24 


UZ. u.zo.uo 


— U ( .ou.zo.ou 


23.64 


22.08 


20.89 


1 nc;-|-n no 
J- . uozcu. uz 


3.7 


21.0 


2.7 




Wl 


SA25 




— 07-91 -m "^9 

U 1 . Z/ J. . U -L . OZ/ 


21.50 


20.06 


19.29 


n 4Q-|-n 09 
u.^ozcu.uz 


2.8 


11.7 


2.0 




Wl 




no.1 fi-nzL fifi 

UZ. lU.U^.UU 


nQ.Qc;.n« ar: 

— Uc^.OO.UU.UO 


91 fiR 


20 26 


1 Q OQ 


n fiQ-Un no 

U.Uc^ZCU.UZ 


16 4 


Ri 

01. u 


9 7 

Z. ( 


V 


Wl 

VV 1 


Q A o'vmG 
oAZ ("'■^ 


no. 1 a. OA no 
Uz: iD:z4.Uo 


nn. ccy. Qn nn 

— uy:o / :oy.uy 


1 7 OQ 

i / .zy 


1 <5 Q c; 
ib.oO 


1 c: Qn 

io.yu 


n 1 Qj_n no 
U. iozLU.Uz 


Q 
Z.O 


c: n 

O.y 


Q Q 
0.0 


A 

A 


W i 


A 00 

dAzo 


no. 1 . Q 1 in 

Uz:lD:ol.iy 


nT. Q 1 . C7 1 Q 

— U7:ol:o7.1o 


21.90 


20.89 


19.80 


n o^^ _Ln n /i 
U.ob±U.U4 


2.4 


12.9 


2.3 


A 

A 


Wl 


c A on 

oAzy 


no. '\ a. Ad OA 
Uz:1d:4d.o4 


nn. 1 . 1 T/i 

— Uy:lo:lb.74 


21.73 


20.53 


19.48 


n TO-Ln nQ 
U.7z±U.Uo 


2.4 


12.1 


2.5 


V 


Wl 


c A on 


no. 1 /( n occ 
Uz: io:4y.zo 


n'y.nQ.OQ on 
— U7:Uo:zo.oU 


21.04 


19.54 


18.85 


n A c^_Ln no 
U.4o±U.Uz 


5.6 


22.6 


Q 

0.7 


A 

A 


W 1 


C A Q1 TTl 


no. 1 '7 . OQ 'y^^ 
Uz:17:zc).7b 


1 n. 1 . ccn Qn 
— lU:lo:oU.oU 


18.68 


17.53 


17.03 


n o'T-Ln no 
U.z7±U.Uz 


3.2 


9.3 


2.3 


^ T 

V 


W 1 


Q A QO 

oAoz 


no. 1 7. on 
Uz: i / :oy. Ob 


1 n. QQ . 1 n nQ 
— iU:oo:iy.yo 


00 Q/l 

zz.y4 


00 nc; 
ZZ. UO 


01 Q/l 

zi.U4 


1 nc;_Ln nc; 
i.UOzLU.UO 


1 n 

i.y 


1 n Q 
iU.o 


Q 
Z.O 


A 

A 


W i 


Q A Q.Q.rn 
oAoo 


Uz. io.u < .zy 


nc . 1 r:.o£i 1 
— UO. iO. ob. ib 


ZZ. 04 


Zi.Zi 


on 1 Q 
zu. io 


n /i O-Un nQ 
U.4ZztU.Uo 


A 
z.4 


Q Q 

y.o 


Q 7 
0. < 


V 


VV i 




no.1 Q.I 4 QQ 

uz. ±0. -L^.OJ:/ 


— lU.UU.UZ.OU 


21 20 


90 ^7 
zu. ( 


90 97 
zu. Z ( 


4fi-l-0 0*^ 
u.^uzcu.uo 


9 R 
Z.O 


11 4 


9 
Z.U 




Wl 

VV 1 


SA35 


uz. J.y.Uc/.OU 


— 04-01 -4*^ "^9 

U^.U-L.^O.OZ/ 


21.43 


19.95 


19.27 


n /I c:_|_n no 
u.^ozcu.uz 


4.3 


17.3 


2.3 




Wl 




no.1 Ot-K(\ 49 
uz. ly.ou.^z 


nc;.o7.c;Q 01 
— uo.z ( .oy.zi 


90 4R 


.00 


1 R R4 


n Qc:_|_n n^ 
u.oozcu.u^ 




O.U 


10 4 




O.U 


V 


Wl 

VV 1 


SA37 


n9-9n-4'^ 1 1 


— 1 0-1^9-1 


22.81 


21.68 


20.48 


n 7q-|-n nQ 

U. ( c/ZCU.UO 


2.2 


11.5 


2.7 




Wl 


SA38 


n9-9n-^fi 4*^ 

uz. zu. ou.^o 


-07-4*^-1 1 71 
— U ( .^0. J. J. . ( J. 


22.91 


21.56 


20.51 


n 71 _|_n nQ 

U. ( J-ZCU.UO 


2.4 


12.1 


2.3 




Wl 


SA39 


UZ.Z J..O J.. J.O 


— 0fi-47-'^9 fifi 

— UU.^( .oz. uu 


21.34 


20.18 


19.16 


n 79-1-0 0*^ 

U. ( ZZCU. UO 


5.2 


26.3 


4.0 


Y 


Wl 


SA40 


n9-9'^-1 41 


— 0fi-9Q-0fi 40 
uu.zy.uu.^u 


21.20 


20.02 


19.21 


^^-\-r\ ofi 

U.OOZCU.UU 


1.9 


8.5 


3.0 




Wl 


SA41 


UZ/.Z/O. -LO.OO 


— 1 O-^R-dR 4R 


21.57 


20.29 


19.53 


p;9-|-0 04 
u.ozzcu.u^ 


6.1 


26.6 


2.0 




Wl 


SA42 


n9-94-nn q9 

UZi.Z(^.UU.c/Z/ 


— 0'^-4R-9'^ S*^ 

UO.^U.ZtJ.OO 


23.12 


22.04 


20.92 


qs-l-0 0^^ 

U.c/OZCU.UO 


2.6 


14.5 


2.7 




Wl 


SA43 


n9-94-n'^ ni 

UZ.Z^.UtJ.UJ- 


— 04-47-07 00 

— u^.^ ( .u ( .uu 


20.00 


18.68 


18.10 


Qfi-UO 04 

U.OUZCU.U^ 


4.3 


15.1 


2.0 




Dl 


SA44 


n9-94-'^4 Qfi 
uz.z^.o^.yu 


— 04-1 1 -"^^ 09 

— u^. J. J. . 00. uz 


24.01 


23.03 


22.06 


fiS-l-0 04 

U.UOZCU.U^ 


1.9 


9.4 


2.7 


Y 


Dl 


SA45 


n9-94-'^'^ 9fi 
uz.z^.oo.zu 


— 04-01 -^7 Sfi 

— U^.UJ-.O ( .ou 


22.49 


21.21 


20.13 


1 1 Q-i-0 07 

J- . J-OZCU.U ( 


3.5 


20.1 


2.0 


Y 


Dl 


SA46 


n9-94-'^Q nfi 

UZ.Z^.Oc/.UU 


— 04-00-4^=; 1 fi 

— U^.UU.^O. -LU 


20.75 


19.33 


18.62 


4Q-I-0 0^=; 

U.^OZCU.UO 


3.2 


12.6 


3.0 


Y 


Dl 


SA47 


n9-94-f^Q 9"^ 

UZ/.Z/^.tJc/.Z/O 


— 04-01 -0*^ 77 

U^.U-L.UO. ( ( 


24.05 


22.83 


21.70 


SO-l-0 04 
u. ouzcu. u^ 


1.9 


10.0 


3.0 


Y 


Dl 


Q A AQzG 


no.oc;. 1 1 n/i 
Uz.zo. ii.U4 


n/1 . K/l . QQ KA 

— U4. 04. 00. 04 


1 Q 70 
io. / Z 


i / .0 / 


1 7 nQ 
i / .Uo 


n QQ-Un ni 
U.ooztU.Ui 


Q 
Z.O 


Q Q 

y.o 


Q 7 
0. / 


A 
A 


T^l 
Ui 


Q A /in 


nO.OK.QQ '7 A 

Uz:zo:oo. / 4 


n/i .nQ.on Qa 
— U4:Uo:zU.ob 


00 nn 

zz.uy 


on RA 
zU.b4 


1 n c:o 
iy.Oz 


U.bzzLU.Ub 


A Q 
4.0 


on A 
zU.4 


n 
z.U 


A 

A 


T^l 
Ui 


c A c^n 


no. OV^. Ad. 1 Q 

Uz:zo:4b. io 


n'T. Q'T. QQ ceo 

— U7:o7:oo.oz 


20.99 


19.50 


18.60 


n c^ /I -Ln no 
U.o4±U.Uz 


5.8 


25.8 


4.0 


^ T 

V 


W 1 




no.o^^.n'y i c^ 
Uz:zb:U7.1o 


n A . OT. o/^ 
— U4:z ^ .Zb.Zb 


19.36 


18.44 


17.97 


n 1 'T-Ln n c^ 
U.17±U.Uo 


Q 

0.7 


c: 

7.0 


2.0 


^ T 
V 


T^l 
Ui 


C A ceo 

dAoz 


no.o'y.on oo 
Uz:z ( :zU.zz 


r\'7 . A n. on 1 n 
— U7:4y:zU.iy 


21.31 


19.81 


19.00 


n ccQ-Ln nQ 
U.Oo±U.Uo 


2.1 


9.3 


2.0 


T T 
V 


W 1 




no. 07. Kn 01 
Uz:z / :oy.zi 


nn.n'y.on q<5 
— uy:U / :zy. ob 


on 71 
zU. / i 


1 n 0^7 

iy.o / 


1 n /1 

iy.4z 


n c: cc-Un nQ 
U.OOzLU.Uo 


Q n 

o.y 


1 7 c; 
i / .0 


n 
z.U 


A 

A 


W i 


c A c: /I 


no. 00 . QO n c 
Uz:z8:oz.U5 


nn. /I n. /I c; /i /i 
— Uy:4y:4o.44 


20.25 


18.74 


18.07 


n A c; -Ln no 
U.4o±U.Uz 


6.3 


25.4 


2.7 


V 


Wl 


Q A c;c; 


no. OQ. 1 7 Q<^ 

Uz:zy:i ^ .ob 


nc;.c:/i.nc: c:/i 
— U0:04:U0.04 


1 Q 

iy.b I 


1 Q Qn 
io.oU 


1 7 'TQ 


n QQ-un nQ 
U.oozLU.Uo 


<^ 
z.b 


Q c: 

y.o 


Q n 
O.U 


A 

A 


W i 


SA56 


02:30:39.96 


-03:50:28.06 


22.36 


21.89 


21.32 


0.27±0.02 


2.1 


6.1 


2.0 


A 


Wl 


SA57 


02:31:06.46 


-05:55:04.63 


21.68 


20.20 


19.46 


0.52±0.03 


3.7 


16.1 


2.0 


A 


Wl 


SA58 


02:32:23.77 


-08:50:38.37 


22.05 


20.76 


20.09 


0.46±0.04 


2.6 


10.6 


2.0 


A 


Wl 


SA59 


02:33:07.05 


-04:38:38.21 


21.56 


20.66 


19.62 


0.79±0.03 


1.9 


9.9 


2.0 


A 


Wl 


SA60 


02:35:01.61 


-09:58:32.76 


21.78 


20.24 


19.07 


0.70±0.03 


4.7 


23.5 


2.7 


A 


Wl 


SA61 


08:48:23.66 


-04:07:15.29 


21.17 


19.63 


18.85 


0.51±0.02 


7.4 


32.0 


2.7 


A 


W2 


SA62 


08:50:07.72 


-01:23:53.30 


22.23 


20.92 


20.35 


0.37±0.04 


3.5 


12.5 


2.7 


A 


W2 


SA63 


08:52:07.18 


-03:43:16.28 


20.91 


19.33 


18.61 


0.48±0.02 


5.0 


20.9 


3.3 


V 


W2 


SA64 


08:52:08.36 


-04:05:28.36 


21.62 


20.26 


19.59 


0.43±0.03 


2.4 


9.4 


2.0 


A 


W2 


SA65 


08:54:25.14 


-03:14:53.11 


24.89 


24.53 


22.84 


0.98±0.10 


1.9 


10.6 


2.7 


A 


W2 


SA66^ 


08:54:46.55 


-01:21:37.08 


19.34 


17.88 


17.28 


0.48±0.11 


4.8 


20.1 


4.0 


A 


W2 


SA67 


08:55:59.92 


-04:09:17.76 


21.06 


19.60 


18.90 


0.45±0.02 


2.1 


8.5 


3.0 


A 


W2 


SA68 


08:57:26.91 


-02:42:26.64 


19.91 


18.44 


17.80 


0.42±0.02 


2.8 


10.8 


2.3 


A 


W2 


SA69 


08:57:35.96 


-01:01:12.55 


21.03 


20.57 


20.29 


0.05±0.23 


2.4 


1.6 


2.3 


A 


W2 



The SARCS sample 



ID Ra Dec g r i ^^phot Rank Type Field 

hms dms mag mag mag h~^Kpc 



SATO 


08: 


:57: 


:49.10 


-01: 


:13: 


:00, 


,73 


19, 


,99 


18, 


,81 


18.26 


0.29+0, 


,03 


3.9 


11.9 


2, 


,0 


A 


W2 


SA71^ 


08: 


:58: 


:48.83 


-02: 


:39: 


:25, 


.79 


19, 


.16 


18, 


.23 


17.83 


0.36+0, 


,10 


3.7 


13.0 


3, 


.0 


V 


W2 


SA72 


08: 


:59: 


:14.55 


-03: 


:45: 


:14, 


,85 


22, 


,01 


20, 


.75 


19.67 


0.74+0, 


,03 


4.5 


23.0 


4, 


.0 


A 


W2 


SA73^ 


08: 


:59: 


:54.54 


-01: 


:32: 


:13, 


.39 


20, 


.87 


19, 


.47 


18.85 


0.66+1, 


.06 


4.3 


21.0 


2, 


.0 


A 


W2 


SA74 


09: 


:00: 


:50.10 


-02: 


:30: 


:54, 


.15 


20, 


.52 


19, 


.23 


18.65 


0.36+0, 


.02 


3.2 


11.3 


2, 


.0 


V 


W2 


SA75 


09: 


:02: 


:20.42 


-02: 


:30: 


:57, 


.28 


22, 


.76 


21, 


.97 


21.41 


0.63+0, 


,04 


3.0 


14.4 


2, 


.3 


A 


W2 


SA76^ 


09: 


:04: 


:07.97 


-00: 


:59: 


:52, 


.85 


22, 


.03 


21, 


.18 


20.22 


0.77+0, 


.04 


2.4 


12.4 


3, 


.3 


A 


W2 


SA77 


09: 


:05: 


:29.75 


-02: 


:03: 


:17, 


.70 


20, 


.94 


19, 


.47 


18.80 


0.42+0, 


.02 


2.4 


9.3 


2, 


.0 


A 


W2 


SA78^ 


09: 


:59: 


:39.17 


+02: 


:30: 


:43. 


.98 


22, 


.72 


21, 


.20 


19.92 


0.74+0, 


.06 


3.2 


16.4 


3, 


.0 


V 


D2 


SA79 


09: 


:59: 


:42.43 


+02: 


:29: 


:56. 


.10 


23, 


.10 


21, 


.64 


20.32 


0.76+0, 


.04 


3.5 


18.1 


3, 


.0 


V 


D2 


SA80 


09: 


:59: 


:55.98 


+02: 


:19: 


:01, 


.79 


23, 


.30 


22, 


.19 


21.06 


1.00+0, 


.04 


2.4 


13.5 


3, 


.3 


A 


D2 


SA81 


10: 


:01: 


:33.74 


+02: 


:21: 


:35. 


.35 


21, 


.96 


20, 


.71 


20.08 


0.69+0, 


.05 


3.0 


14.9 


2, 


.0 


V 


D2 


SA82 


10: 


:01: 


:47.79 


+02: 


:22: 


:06. 


.55 


22, 


.92 


21, 


.44 


20.20 


0.69+0, 


.05 


3.5 


17.4 


2, 


.7 


V 


D2 


SA83^ 


10: 


:02: 


:11.22 


+02: 


:11: 


:39. 


.46 


23, 


.64 


22, 


.02 


20.77 


0.89+0, 


.05 


2.6 


14.1 


4, 


.0 


V 


D2 


SA84 


10: 


:02: 


:11.67 


+02: 


:29: 


:55. 


.24 


22, 


.81 


21, 


.55 


20.56 


0.77+0, 


.05 


1.9 


9.9 


3, 


.0 


A 


D2 


SA85 


10: 


:02: 


:14.85 


+02: 


:37: 


:36, 


.47 


23, 


.51 


22, 


.20 


21.15 


0.65+0, 


.05 


2.1 


10.2 


3, 


.0 


V 


D2 


SA86 


13: 


:56: 


:49.33 


+55: 


:27: 


:07. 


.00 


20, 


.41 


18, 


.86 


18.21 


0.46+0, 


.03 


3.7 


15.1 


2, 


.5 


A 


W3 


SA87 


13: 


:57: 


:25.48 


+53: 


:17: 


:43, 


.96 


20, 


.52 


19, 


.03 


18.15 


0.54+0, 


.02 


3.5 


15.6 


3, 


.0 


A 


W3 


SA88 


13: 


:59: 


:47.26 


+55: 


:35: 


:37, 


.57 


23, 


.40 


22, 


.63 


21.57 


0.87+0, 


.04 


2.2 


11.9 


2, 


.3 


A 


W3 


SA89 


14: 


:00: 


:40.17 


+56: 


:07: 


:49, 


.41 


20, 


.47 


19, 


.08 


18.48 


0.42+0, 


.03 


3.7 


14.3 


3, 


.0 


A 


W3 


SA90 


14: 


:01: 


:10.46 


+56: 


:54: 


:20, 


.51 


20, 


.89 


19, 


.42 


18.57 


0.53+0, 


.03 


3.7 


16.3 


3, 


.7 


A 


W3 


SA91 


14: 


:01: 


: 44.90 


+53: 


:02: 


:09, 


.62 


22, 


.01 


20, 


.50 


19.61 


0.56+0, 


.03 


3.0 


13.6 


3, 


.7 


A 


W3 


SA92^ 


14: 


:01: 


:56.39 


+55: 


:44: 


:46. 


.78 


20, 


.79 


19, 


.40 


18.70 


0.50+0, 


.03 


2.8 


12.0 


2, 


.0 


V 


W3 


SA93 


14: 


:02: 


:47.90 


+57: 


:08: 


:52, 


.04 


23, 


.70 


22, 


.75 


21.68 


1.22+0, 


.04 


3.2 


18.6 


2, 


.0 


A 


W3 


SA94 


14: 


:03: 


:51.68 


+57: 


:23: 


:50, 


.41 


24, 


.08 


22, 


.48 


21.72 


0.51+0, 


.03 


0.0 


0.0 


2, 


.7 


A 


W3 


SA95^ 


14: 


:04: 


:54.46 


+52: 


:00: 


:24, 


.70 


20, 


.10 


18, 


.66 


17.88 


0.49+0, 


.03 


2.2 


9.3 


2, 


.7 


A 


W3 


SA96^ 


14: 


:05: 


:54.33 


+54: 


:45: 


:48, 


.68 


19, 


.50 


18, 


.18 


17.49 


0.41+0, 


.03 


2.8 


10.7 


3, 


.0 


V 


W3 


SA97 


14: 


:08: 


:13.82 


+54: 


:29: 


:08, 


.12 


20, 


.28 


18, 


.79 


18.04 


0.48+0, 


.02 


8.0 


33.4 


4, 


.0 


A 


W3 


SA98 


14: 


:11: 


:20.53 


+52: 


:12: 


:09, 


.91 


20, 


.16 


18, 


.75 


17.93 


0.52+0, 


.03 


18.4 


80.3 


2, 


.0 


A 


W3 


SA99^ 


14: 


:13: 


:55.43 


+53: 


:43: 


:44. 


.72 


19, 


.03 


17, 


.83 


17.26 


0.29+0, 


.03 


2.4 


7.3 


2, 


.0 


A 


W3 


SAIOO 


14: 


:14: 


:47.19 


+54: 


:47: 


:03, 


.59 


21, 


.19 


19, 


.67 


18.45 


0.63+0, 


.02 


14.7 


70.3 


3, 


.7 


A 


W3 


SAlOl 


14: 


:16: 


: 44.52 


+56: 


:42: 


:16. 


.18 


22, 


.99 


21, 


.41 


19.94 


1.29+0, 


.16 


3.5 


20.5 


2, 


.0 


A 


W3 


SA102^ 


14: 


:19: 


:12.17 


+53: 


:26: 


:11. 


.44 


21, 


.83 


20, 


.30 


19.11 


0.69+0, 


.02 


9.9 


49.2 


3, 


,7 


A 


W3 


SAIOS^^ 


14: 


:19: 


:17.25 




:17: 


:28. 


.63 


20, 


.78 


19, 


.50 


18.72 


0.47+0, 


.03 


4.1 


16.9 


3, 


.0 


V 


W3 


SA104 


14: 


:21: 


:02.56 


+52: 


:29: 


:42, 


.51 


17, 


.74 


16, 


.79 


16.33 


0.18+0.1 


01^ 


11.7 


24.9 


2, 


.0 


A 


D3 


SA105 


14: 


:21: 


:18.35 


+52: 


:50: 


:22, 


.37 


21, 


.40 


19, 


.89 


19.14 


0.47+0, 


.05 


2.8 


11.6 


3, 


.0 


V 


D3 


SA106 


14: 


:22: 


:58.34 


+51 ■ 


:24: 


:39. 


.50 


22, 


.78 


21, 


.75 


20.80 


0.74+0, 


.04 


1.9 


9.7 


2, 


.5 


A 


W3 


SA107 


14: 


:23: 


:49.27 


+57: 


:26: 


:33. 


.90 


23, 


.20 


22, 


.15 


21.21 


0.69+0, 


.06 


2.2 


10.9 


2, 


.0 


A 


W3 


SA108 


14: 


:25: 


: 44.27 


+57: 


:07: 


:24, 


.47 


25, 


.45 


23, 


.87 


22.53 


0.86+0, 


.04 


4.5 


24.2 


2, 


.7 


A 


W3 


SA109 


14: 


:26: 


:08.04 


+57: 


:45: 


:23. 


.90 


20, 


.56 


19, 


.49 


18.99 


0.39+0, 


.03 


3.2 


11.9 


2, 


.0 


A 


W3 


SAllO^ 


14: 


:28: 


:10.54 


+56: 


:39: 


:48. 


.36 


17, 


.67 


16, 


.76 


16.30 


0.80+0, 


.28 


4.1 


21.5 


2, 


.3 


A 


W3 


SAlll 


14: 


:28: 


:34.82 


+52: 


:13: 


:06, 


.44 


22, 


.28 


20, 


.75 


19.94 


0.52+0, 


.03 


5.0 


21.8 


2, 


.7 


A 


W3 


SA112 


14: 


:30: 


:00.65 


+55: 


:46: 


:47, 


.97 


21, 


.58 


20, 


.02 


19.12 


0.55+0, 


.02 


4.3 


19.3 


4, 


.0 


A 


W3 


SA113 


14: 


:31: 


:39.77 


+55: 


:33: 


:22, 


.81 


21, 


.83 


20, 


.57 


19.42 


0.71+0, 


.03 


3.0 


15.1 


3, 


.0 


V 


W3 


SA114 


14: 


:31: 


:52.67 


+57: 


:28: 


:36, 


.73 


22, 


.87 


21, 


.52 


20.19 


0.83+0, 


.03 


3.5 


18.6 


2, 


.3 


A 


W3 


SA115 


14: 


:34: 


:03.87 


+51: 


:21: 


:36, 


.07 


19, 


.98 


18, 


.66 


18.06 


0.39+0, 


.02 


2.6 


9.6 


2, 


.0 


A 


W3 


SA116 


14: 


:34: 


:34.69 


+56: 


:59: 


:20, 


.17 


21, 


.19 


19, 


.63 


18.74 


0.57+0, 


.02 


4.1 


18.7 


2, 


.7 


A 


W3 


SA117^ 


22: 


:01: 


:51.79 


+04: 


:10: 


:08, 


.42 


18, 


.53 


17, 


.73 


17.19 


0.43+0, 


.04 


7.3 


28.7 


2, 


.7 


A 


W4 


SA118^ 


22: 


:02: 


:01.66 


+01: 


:47: 


:09, 


.57 


18, 


.82 


17, 


.63 


17.07 


0.30+0, 


.02 


5.0 


15.6 


3, 


.0 


A 


W4 


SA119^ 


22: 


:03: 


:29.03 


+02: 


:05: 


:18, 


.89 


21, 


.24 


19, 


.99 


19.37 


0.38+0, 


.04 


2.6 


9.5 


4, 


.0 


V 


W4 


SA120^ 


22: 


:05: 


:06.92 


+01: 


:47: 


:03. 


.71 


21, 


.20 


19, 


.91 


19.15 


0.46+0, 


.06 


2.1 


8.6 


2, 


.0 


A 


W4 


SA121 


22: 


:06: 


:42.03 


+04: 


:11: 


:30, 


.85 


21, 


.20 


19, 


.81 


18.88 


0.62+0, 


.03 


3.7 


17.6 


3, 


.0 


A 


W4 


SA122 


22: 


:13: 


:06.93 


-00: 


:30: 


:37. 


.05 


21, 


.19 


19, 


.98 


18.81 


0.69+0, 


.02 


2.8 


13.9 


3, 


.0 


A 


W4 


SA123 


22: 


:13: 


:31.85 


+00: 


:48: 


:36, 


.14 


23, 


.37 


21, 


.87 


20.56 


1.00+0, 


.03 


4.8 


26.9 


4, 


.0 


A 


W4 


SA124 


22: 


:14: 


:09.57 


-17: 


:30: 


:56, 


.23 


22, 


.63 


21, 


.08 


19.85 


0.83+0, 


.05 


7.4 


39.4 


4, 


.0 


V 


D4 


SA125 


22: 


:14: 


:18.82 


+01: 


:10: 


:33, 


.85 


20, 


.31 


19, 


.24 


18.84 


0.74+0, 


.08 


0.0 


0.0 


3, 


.0 


A 


W4 


SA126 


22: 


:17: 


:29.38 


-00: 


:38: 


:36. 


.60 


19, 


.93 


18, 


.94 


18.32 


0.78+0, 


.02 


1.9 


9.9 


2, 


.0 


A 


W4 


SA127 


22: 


:21: 


:43.74 


-00: 


:53: 


:02, 


.89 


19, 


.36 


17, 


.98 


17.35 


0.39+0, 


.02 


4.7 


17.4 


3, 


.3 


V 


W4 



m: This galaxy falls within the masked region as per the catalog from which the magnitudes and the redshift are 
extracted, z: The magnitudes and/or redshift are not from the Coupon et al. catalog instead are measured by 
the author using sextractor and/or zebra (Feldmann et al. 2006), respectively. C: This lens is identified in both 
D2 and COSMOS fields. Note that other lenses within D2 have not been reported in the COSMOS lens sample 
(Faure et al.. .2008. V R: This lens is also found in the RCS (see G03). G: This lens is also a part of the SL2S-RINGS 
sample. 



3. RESULTS AND DISCUSSION 

In the following subsections, we describe the main find- 
ings from the SARCS sample and constraints on average 
properties of the lens population using statistical prop- 
erties of the arcs. Note that many of the lens candidates 
are not confirmed lenses yet and hence, the results should 
be taken as indicative. Firstly, based on the photometric 



redshift s of the lenses, we study the lens redshift distri- 
bution. Subsequently, we discuss about the abundance 
of giant arcs and presence of radial arcs in the sample. 
Finally, we measure the azimuthal distribution and im- 
age separation distribution of the arcs and argue their 
importance as diagnostics for understanding the matter 
distribution of the lenses, statistically. 
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6 020320-073451 I 7 020350-094254 I 9 020503-110547 I 10* 020649-(«5701 



14 020929-06431 




21 021325-074355 I 22* 021409-053532 
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23 02141 K040503 I 27 021624-095739 1 30 021649-070324 021807-051536 I 36 02?956-052^59 





46 022439-040045 




48 022511-045434 




55 022917-0554061 63 085207-034^161 66 085447-0121371 67 085600-040918 



71 085849-023926 




76 O904O8-OO59531 78 O95939+O230441 79 "^95942+022956 



Fig. 3. — The SARCS sample showing the 54 promising lens systems with rank of 3 and above. All color cutouts are made from CFHT 
imaging in r and i bands. The cutouts are ~30^^on the side. The candidate SA78 is known as COSMOS 5939+3044. 
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123 221332+004836 1124 221440-17*056 
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127 222144-0t)5303 



Continue Fig. [3 The candidates SA83 and SA102 are known as COSMOS 0211+1139 and RCS 1419.2+5326, respectively. 
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Fig. 5. HST /-band (F814W) imaging of the SARCS lens candidates within the D2 field extracted from the COSMOS archive. The 
IDs and co-ordinates are labelled for each system. The arrows point towards the putative lensed arcs. The arc in SA79 is too faint to 
be visible in the /-band HST image shown here. All candidates are new detections except SA78 (COSMOS ID - 5939+3044) and SA83 
(COSMOS ID - 0211+1139). The images are 6^^on the side. North is up and East is left. 



TABLE 3 

SARCS CANDIDATES WITH FOLLOW-UP INFORMATION FROM THE SL2S COLLABORATION. 



RA 


DEC 


Reference 




Comment 


02:09:57.67 


-03:54:57.08 


VM 




V, 


H 


02:13:24.52 


-07:43:54.82 




0.72, - 


K, 


V, H 


02:14:08.07 


-05:35:32.39 


L09,V11 


0.444, 


V, 


H 






VM 


1.023+0.001/1.7+0.1 






02:14:11.24 


-04:05:02.71 


M11,VM 


0.609, - 


K, 


V, H 


02:16:49.25 


-07:03:23.80 


VM 




V, 


H 


02:18:07.29 


-05:15:36.16 


Mil 


0.647, - 


V, 


H 


02:19:56.42 


-05:27:59.21 


VM 




V, 


H 


02:21:51.18 


-06:47:32.66 


L09,VM 


0.618, 


V, 


H 


02:25:11.04 


-04:54:33.54 


Rll 


0.238, 1.199 


K, 


H 


02:25:46.13 


-07:37:38.52 


L09 


0.511, 


G 




08:52:07.18 


-03:43:16.28 


VM 




V, 


H 


08:54:46.55 


-01:21:37.08 


L09,L10 


0.3530+0.0005,1.2680+0.0003 


K, 


V, H 


08:58:48.83 


-02:39:25.79 






H 




08:59:14.55 


-03:45:14.85 


L09,VM 


0.647, - 


V, 


H 


09:04:07.97 


-00:59:52.85 






H 




10:02:11.22 


+02:11:39.46 






H 




14:08:13.82 


+54:29:08.12 


L09 


0.416, - 


s, 


H 


14:14:47.19 


+54:47:03.59 






H 




14:19:12.17 


+53:26:11.44 






H 




14:19:17.25 


+51:17:28.63 






H 




14:30:00.65 


+55:46:47.97 


TIO 


0.497+0.001, 1.435+0.001 


G, 


H 


14:31:39.77 


+55:33:22.81 


TIO 


0.669+0.001, - 


G, 


H 


22:03:29.03 


+02:05:18.89 






H 




22:13:06.93 


-00:30:37.05 


L09 


■ 


H 




22:13:31.85 


+00:48:36.14 


L09 




H 




22:14:18.82 


+01:10:33.85 






H 




22:21:43.74 


-00:53:02.89 


L09 


0.334, - 


s, 


H 



Note. — Col 3: L09- Limousin et al." ("2009), LIO- Limousin et al." fMTO), TlOJThaniavur et al." ("20101), Rll-"Rufi"et al." (l20Tll), Vll- 
IVerdugo et al. (2011) and Mll-Muhoz et al. (2011, in prep, PI- V. Motta, ESO-080.A-0610). VM is assigned to those systems which are 
being followed up for spectroscopy with the VLT by V. Motta (PI, ESO-086.A-0412). Col 4: Spectroscopic redshifts for the lens galaxy 
and lensed source from the follow-up observations along with their error bars, if available. Col 5: V-VLT, K-Keck or G-Gemini indicates 
telescopes used for follow-up spectroscopy and S-SDSS spectroscopy. H corresponds to HST imaging followed up by J-P. Kneib (PI, C15 
and C16) and/or R. Gavazzi (PI, C17). 
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Fig. 6. — Redshift distributions of the lens galaxies from the 
SARCS sample. The solid histogram corresponds to lenses from 
the whole sample whereas the dashed histogram corresponds to 
lenses with giant arcs only. The respective means of their redshift 
distributions are z = 0.58 ± 0.22 and z = 0.64 ± 0.19 where the 
uncertainties indicate 68% confidence level. 
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3.1. Lens redshift distribution 

The lenses producing giant arcs (e.g. length-to-width 
ratio > 10) consist of clusters which lie at the high end 
of the halo mass function. The redshift distribution and 
abundance of such massive clusters depend on the param- 
eters of a cosmological model. Bartelmann et al. (1998) 
estimated that most of the lensing clusters giving rise to 
giant arcs should peak at around z ^ 0.3 — 0.4 for the 
currently accepted standard cosmological model. 

We use the CFH TLS photometric redshift catalogs 
(jCouDon et al.ll2009D wh ich are generated from the code 
LE PHARE filbert et al.l l2006[ ) to calculate the SARCS 
lens redshift distribution. The accuracy on the redshifts 
of galaxies in the WIDE with magnitudes i < 21.5 is 
^Az/ii-\-z) ^ 0.037 and with magnitudes 22.5 < i < 23.5 
is (Jaz/{i-\-z) ^ 0.08. In Fig.[6l we show the redshift dis- 
tribution for all the lenses in the SARCS sample (solid 
histogram) and for lenses consisting of giant arcs only 
(dashed histogram). In the appendix B, we describe how 
we estimated the mean and la uncertainty given the red- 
shift measurement errors (see Table [2]). We find that the 
mean of the lens redshift distribution for the SARCS 
sample is z = 0.58 ± 0.22 and that for the sample of gi- 
ant arcs is z = 0.64 ±0.19. We note that the mean of 
the giant arcs sample is a t a higher redshift compared 
to the peak expected from iBartelmann et al.l ([l998) but 
certainly consistent within 2a confidence interval. For 
comparison, the RCS sample G03 finds that most of the 
lenses with giant arcs have redshifts upwards of ~ 0.6. 

We make a qualitative comparison of the redshift dis- 
tributions of the lens populations from other surveys in 
the literature (not shown in the figure). However, we 
note that these surveys have significantly different selec- 
tion functions and hence, quantitative conclusions should 
not be drawn. The distribution of lens sample from C AS- 



TLE^3 peaks between 0.3 and 0.4 which is lower than 
the mean redshift of the SARCS sample. CASTLES 
has a large enough sample of lenses but is not homo- 
geneously selected. Nevertheless, the peaks are consis- 
tent within 2a (assuming an error of 0.05) in spite of the 
differenc es in the sample selection. The COSMOS lens 
sample (|Faure et al.ll2008[ ), on the other hand, is fairly 
homogeneous but the sample size of confirmed lenses is 
relatively small. COSMOS has a bimodal lens redshift 
distribution with a minimum at ~ 0.5 — 0.6. This is in 
stark contrast with the redshift distribution of SARCS 
(see Fig. [6]). 

3.2. Giant and radial arcs 

Here, we report detections of giant and radial arcs from 
our sample. The arcs, both radial and tangential, pro- 
duced in massive clusters are conventionally referred to 
as giant arcs, if their length-to- width (l/w) ratio is larger 
than about 8 or 10. For a (non-singular) circularly sym- 
metric projected density profile, the inverse magnifica- 
tion matrix of a lensed image has two eigenvectors, one 
in the radial and the other in the tangential direction. If 
the radial eigenvalue becomes 0, then the arcs are mag- 
nified and stretched radially with respect to the center 
of the lens and are called radial arcs. If the tangential 
eigenvalue becomes 0, then the arcs are magnified and 
stretched tangentially and are called tangential arcs. 

3.2.1. Giant arcs 

The statistics of giant arcs allow detection of clus- 
ters at the massive end of halo mass function that 
are rare. The statistics of such rare massive struc- 
tures is sensitive to the cosmological model of the Uni- 
verse. Hence, several attempts have been made to predict 
the giant arcs statistics (e.g., B artelmann et al. 1998; 
IWambsganss et al.l 120041 : iDalal et al.l 12004 henceforth, 
D04) that suggested a large discrepancy compared to 
the observed abund ance of arcs from a well-defined clus- 
ter p opulation fe.^..|L uDDino et al."1999': 'Gonzal ez et al.l 
\mm iGladders er^^003; Li et al. 2005, 2006). How- 
ever, the discrepancy has been substantially diminished 
due to more realistic assumptions such as using a realis- 
tic source redshift distribution and improved predictions 
which consider the contribution of central galaxy or sub- 
structure within the hal os (e.g.. iMeneghetti et al.|[2QQQl 
[200l iHoresh et al.l[2Q05l ). 

We present the abundance of giant arcs in our sam- 
ple which could be tested with predictions from realistic 
simulations by taking into account the observational lim- 
itations. Within a total unmasked area of '^150 deg^, 8 
of the arcs have l/w of ~10 or above. Additionally, 4 
more arcs have l/w ~ 8 which are included in the sam- 
ple of giant arcs because these appear to be broken either 
due to noise or due to a satellite galaxy. Since one of the 
giant arcs is from the DEEP, we use rest of the 11 gi- 
ant arcs from the WIDE data for comparison with RCS. 
We use the primary sample of RCS which has roughly 
similar depth compared to WIDE data. The primary 
RCS lens sample has 6 arcs with l/w > S found within 
a total area of 90 deg^ G03. Therefore, the RCS sample 
has 0.07 ± 0.03 arcs deg~^ which is consistent with the 
0.07 ± 0.02 arcs deg-^ from the SARCS sample. 

■^^ http:/ /www. cfa.harvard.edu/castles/ 
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Fig. 7. The CFHTLS gri image of SL2SJ141447+544704 (SAlOO) on the left and HST F606W {V band) image on the right. The 
high resolution of HST suggests two radial arc candidates labeled as Rl and R2 whereas the tangential arc labeled as Tl is a confirmed 
lensed arc. The CFHTLS image is 74.4^^on the side and the HST image is ~28.3^^on the side. 



3.2.2. Radial arcs 

The radial arcs are formed when the source fahs on the 
radial caustic (see Appendix A). The size of the radial 
caustic and hence, the cross-section to form radial arcs 
depends on the slope of in ner regio ns of the density pro- 
file of the lens (e.g., Hattori et al.L l999). The statistics 
of radial arcs can thus, be used to constrain the slope 
of the central density profiles of clusters and thereby, al- 
low a better understanding of the nature of dark matter 
(iMolikaw a fc Hattoril l2QQlf ) . More thorough investiga- 
tions have been carried out to test the effects of realistic 
assumptions of properties such as lens ellipticity, source 
size and ellipticity on the statistics of rad i al and tan- 
gential arcs (e.g., Keeton 2001; Qgurl l2002f ). iSand et al.l 
1^2005) used archival HST (WFPC2) data on various lens- 
ing cluster samples to do a systematic study of the num- 
ber ratio of radial to tangential arc statistics. They ac- 
counted for the effects due to the central galaxy in the 
expected number ratio and placed loose constraints on 
the slope of the average inner density profile of the dark 
matter. They underlined the need of larger observational 
datasets to further explore effects of substructure, mass 
of the central galaxy and homogeneity of the sample. 

We report detections of radial arc candidates in the 
SARCS sample during the visual inspection. We found 
1 candidate which appears like a radial arc in the system 
SL2SJ141447+544704 (SAlOO in Table E]). This lens sys- 
tem also has a giant tangential arc (Tl) with the same 
color as the radial arc (Rl, see the left panel of Fig. [7]). 
The subsequent HST observation in the F606W band 
shows another radial arc candidate (R2, see right panel 
of Fig.[7|). The nature of the radial arc candidates needs 
to be verified with spectroscopy and mass modeling since 
these could be blue edge-on disk galaxies. Radial arcs are 
usually faint and are overshadowed by the bright central 
galaxies near which they are formed. Nevertheless, we 
argue that a follow-up imaging of the promising SARCS 
candidates at high resolution will help in creating a ho- 
mogeneous sample of radial and tangential arcs which 
could provide crucial insights in understanding the den- 



sity distributions in the inner regions of the lens systems. 
3.3. Galaxy-arcs orientation 

In this section, we quantify the angular distribution of 
arcs with respect to the major axis of the ellipticity of 
the lensing halo. D04 suggested that triaxial dark mat- 
ter h alos, when acting as lenses, usually form lip caustics 
(e.g. JHattori et al]ll999l ) and the formation of giant arcs 
tends to be at the ends of the lip caustics. Using nu- 
merical simulations, they showed that the giant arcs are 
oriented very close to the major axis of the dark matter 
halo. If a central galaxy is further added to the halo, 
then it appeared to isotropize the angular distribution of 
arcs to a small extent. In order to compare these pre- 
dictions to the observations, D04 measured the angular 
distribution of giant arcs from the EMSS cluster sample 
under the assumption that the ellipticity of the lensing 
galaxy is a good representation of the ellipticity of the 
underlying DM halo. They found that most of the giant 
arcs had an orientation of < 45 deg consistent with their 
predictions. 

We use the same definition as that used in D04 for 
the orientation of the arcs, that is, the angle between 
the major axis of the lens galaxy and the line connect- 
ing the center of the lens galaxy to the midpoint of the 
arc. The measurement of orientation of real arcs tends 
to be somewhat subjective since there exists an ambi- 
guity in defining the extent of the arc which is required 
to find the midpoint of the arc. The midpoint of the 
arc and its orientation from the major axis of the lens 
is calculated manually. When there is a single domi- 
nant lensing galaxy we measure the position angle (PA) 
of the major axis of the lens galaxy with sextractor 
(jBertin fc Arnoutslll996f ) otherwise we adopt the follow- 
ing strategy. If there are two or three (nearly collinear) 
lens galaxies with similar colors and brightness, then the 
PA of the ellipticity is assumed to be along the line join- 
ing the lens galaxies. If a circle going through the arc en- 
closes multiple lens galaxies with comparable brightness 
but no obvious elongation, then such system is rejected. 
If one of the lens galaxies is predominantly brighter than 
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the companion galaxies, then the ellipticity of the domi- 
nant galaxy is assumed to reflect the overall ellipticity of 
the gravitational potential and hence, of the caustics. 

After applying the above selection cuts, we are left 
with 36 all arc candidates and 11 giant arcs with an av- 
erage ranking > 2.5. The angular distribution of the 
giant arcs from the SARCS sample is shown in the left 
panel of Fig. [H We compare the observed distribution 
to the expected distributions from D04 which correspond 
to two cases: a) for dark matter only lensing halo (solid 
curve) and b) a dark matter halo to which a galaxy of 
3 X 10^^ h~^MQ is added (dashed curve). The distribu- 
tion of giant arcs from our sample appears to follow the 
expected trend although it is not possible to distinguish 
between the two cases. We also show the distribution of 
orientation of all arc candidates (see left panel of Fig. [8]) 
which seems to follow the expected anisotropy. This ex- 
tends the result obtained by D04 to groups-scale lens 
systems. Furthermore, we split the all arcs sample into 
a small Ra (< 5'') and large Ra (>= 5'') samples to see 
the dependence of the angular distribution on the Ein- 
stein radius (or the arc radius). The choice of cut at 
5''is arbitrary. The baryonic matter tends to make the 
matter distribution in halos more spherical at the cen- 
ter. Therefore, the angular distribution of arcs for the 
sample with smaller Ra is expected to be more isotropic 
compared to the sample with large Ra- However, we do 
not find any clear differences in the angular distribution 
of small and large Ra for the SARCS sample within the 
uncertainties (see the right panel of Fig. [8]). We note 
that these measurements may have some systematic er- 
rors due to the orientation dependence of the selection 
function. Ambiguities also exist in the definition of the 
orientation in few cases, especially when multiple candi- 
date lens galaxies are involved. In addition, the current 
analysis consists of lens candidates as opposed to con- 
firmed lens systems. Therefore, our conclusions should 
be treated more of a qualitative nature. 

3.4. Image separation distribution 

The image separation distribution (ISD) is sensitive 
to the halo mass, abundance of the lens population, 
the mass distribution in the lens and the source red- 
shift. Therefore, the ISD measured from galaxy to clus- 
ter scales contains information about the cosmological 
parameters and various scaling relations between galaxy 
properties and halo mass. Hitherto, the ISD has been 
measured either at small image separations (0) primarily, 
with lens samples such as CASTLE S fe.p^.. iKeeton et al. 
[2000; Khare 2001) and CLASS fe.g. JKochanek &: White 
[200 Ic iOgurii i2006) or at large with cl uster-scale lenses, 
for example, the MACS sample (IZitrin et al. 2011). With 
the SARCS sample, we can probe the intermediate mass 
regime corresponding to group scale lenses. The SARCS 
sample is selected based on the presence of elongated arc- 
like features and is not directly biased towards selecting 
a specific lens population. The ISD measured from the 
lens samples is referred to as the observed ISD and the 
ISD calculated from various models is referred to as the 
expected/predicted ISD throughout this paper. 

3.4.1. Model for the expected distribution 

In the past, the expected ISD was computed with 
models consisting of either galaxy dominated (e.g.. 



iTurner et al .11 19841) or dark matter dominated lensing ha- 
los fe.g.. iNaravan White! Il988f ). However, since the 
density profile of dark matter is known to be affected 
due to the presence of baryons at the center, the need 
for a more complex model to expla i n the observed ISD 
was pointed out by iKeeton et al.l (119981) and demon- 
strate d by subsequent stu dies fe.g., iPorciani fc Madaul 
1200(1 iKochanek White! I200TI : iQguril l2006l henceforth, 
O06). We follow the framework developed in O06 and 
lOguri etld] ()2002f ) to calculate the expected ISD in order 
to compare it with the observed ISD from the SARCS 
sample. We also adopt the cosmology used in O06 
which constitutes of the following cosmological param- 
eters: Qrn = 0.3, Qa = 0.7 and <J8 = 0.9. We describe in 
detail how each of the model parameters are calculated 
before proceeding to the comparison of various models. 

The probability for a source at redshift, Zs to get lensed 
with image separation greater than is given by 

P{>0',zs)= [ ^ dzi^ [ dMn{M,zi) 
Jo CLZ/ Jo 

xaiens e(M-M), (4) 

where x is the comoving distance and zi is the redshift 
of the lens, M corresponds to the halo mass, n{M,zi) 
is the halo mass function, M is the minimum halo mass 
that causes an image separation equal to ^, and 6 is the 
Heaviside step function. The biased lensing cross section, 
(Jiens is measured in comoving units and is given by 

CTlens = 7ry'^^QB{zs) , (5) 

where yr(,o corresponds to the radius of the outermost 
caustic in the lens plan43 and it depends on the mat- 
ter density distribution around the lens. The quan- 
tity B{zs) denotes the magnification bias which causes 
sources, fainter than the limiting magnitude of the sur- 
vey, to be detected in the sample. It is the ratio of the 
number of sources that can be potentially lensed into an 
image with luminosity L to the number of sources that 
have an intrinsic luminosity L. In general, the magnifica- 
tion bias depends upon the luminosity of the source and 
under the assumption of spherical symmetry, it is given 

by 

1 f^'^ / L \ dL 

^(^^'^) = ^^7 7TTF <^2/27ry$ 2:,, -—)—-, 

7Ty^<i>{Zs,L)dL Jo V i^iyj^ n{y) 

(6) 

where ^(^s, L) is the true source luminosity function and 
fi{y) is the lensing magnification at an angular position 
y inside the caustic. Under the assumption of a power- 
law luminosity function which does not evolve with red- 
shift, the luminosity dependence of the magnification 
bias drops out. 

Differentiating both sides of Eq. HI we obtain the dif- 
ferential probability 

dP f^' dy dM 

= rdzz^n(M[^],zOaiens^7,) 

For lens models with spherical symmetry that are considered 
here, the radial caustic is the outermost caustic and a source lying 
within the outermost caustic gets strongly lensed, i.e., multiply 
imaged. 
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Fig. 8. — Angular distribution of arcs with respect to the lens galaxy. The data points in both the panels are from the SARCS sample. 
In the right panel, the "all arcs" sample is divided into small and large Ra- The expected curves for dark matter (DM) only (solid) and 
DM+3 X lO^^/i~^M0 (dashed) are taken from D04 which are calculated from the cluster sample of GIF simulation ( Kauffmann et aOl999I V 



which can be directly related to the ISD of the observed 
lens sample. 

In O06, the above equation is used to predict the ISD 
resulting from different components in a given halo. To 
use the above equation, we have to specify the distribu- 
tion of mass inside a halo and the halo mass function. 
The former allows us to calculate M{0) and the biased 
cross section for lensing. For the latter, O06 assumed 
the then state-of-the-art calibration of the mass function 
given by Sheth fc TormenI (jl999[ ) . For the mass distribu- 
tion inside a halo, O06 considered the following different 
components. 

• At the center of the dark matter halo, the mat- 
ter density is dominated by the central galaxy 
and the total matter distribution is very close 
to that of a singular isothermal sphere (SIS). 
This distribution is given in terms of the one- 
dimensional velocity dispersion (avei) such that 
p(r) (X cr^i/ r^. Following Sect 2.3.1 of 
lOguri et al. (2002), we calculate the size of the ra- 
dial caustic {yr£,o) in comoving coordinates, the to- 
tal magnification of the lensed images and the re- 
lation between the velocity dispersion and the im- 
age separation for the SIS mass distribution. To 
relate the mass of the halo to the velocity dis- 
persion of the central galaxy, O06 used galaxy 
scaling relations. They adopte d the halo ni ass- 
luminosity relation from i Vale fc O striker! ()2004f ) us- 
ing the abundance matching tech nique and the 
Faber - Jackson relation obtained by 'Bern ardi et al.l 
()2003l ) using SDSS galaxies to relate the luminos- 
ity to the veloc ity dispersion of the galaxy. Fol- 
lowing ^ Cooray fc Milosavljevic (2005), they also 
included a log-normal scatter in the halo mass- 
lum inosity relation with a sca tter of 0.25 dex (see 
also [More et al.l[2009bl l2011b[ ). 



• On large scales, the distribution of dark mat- 
ter follows the Navarro- Frenk- White (NFW, 



iNavarro et al.lflQQTl ) profile given by 
p{r) 



{r/rs){l^r/rs)^ ' 



(8) 



where the scale radius = i?vir/c, the concentra- 
tion parameter c is related to the mass M with a 
considerable scatter. O06 used the mean relation 
betwe en concentration and mass of iBullock et al.l 
(l200Th . given by 



10 



1 + z 



M 



-0.13 



(9) 



The distribution of concentrations at fixed halo 
mass is assumed to be log-normal with a scatter of 
0.3 dex. For the NFW case, we calculate the caus- 
tic size, the magnification and the relation between 
halo mass and the image separation numerically. 
We defer the details of our procedure to calculate 
these quantities to the appendix A. 

• In addition to the above two simple profiles, O06 
also considered a combined total profile which in- 
cluded the central galaxy, the dark matter halo and 
the effect of adiabatic contraction (AC) of dark 
matter in response to the baryonic component of 
the galaxy at the center. In this case, O06 as- 
sum ed the central ga laxy to have a Hernquist pro- 
file (|Hernquistlll990[ ). given by 



p{r) 



27r {r/n){n + r)^ 



(10) 



where M5 is the stellar mass of the galaxy and 
Tb is a core radius. The stellar mass M5 was 
obtained using the ha lo mass-luminosity relation 
found by >Vale fc OsFriker (2004) and adopting a 
constant mass-to-light rat io of S.OhjoMG^/ Lc:). The 
scaling relations of B ernardi et al.l ()2003f ) can be 
used to obtain the effective radius Rq as a function 
of the luminosity of the galaxy which is related 
to the core radius of the Hernquist profile such 
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Fig. 9. — Image separation distribution. Theoretically calcu- 
lated image separation distribution curves for SIS profile (dotted), 
NFW profile (dashed) and total profile (solid) following the O06 
model. The data points from the SLACS (green), SARCS sample 
(magenta) and MACS (cyan) where the vertical bars indicate Pois- 
sonian errors and horizontal bars show the bin width. As discussed 
in Sect. 13.4.21 we multiply the ISDs of the lens samples by their 
respective P(> ^cut)- 



that = 0.551i?o- The AC is carr ied out using 
the an alytical formalism presented bv lGnedin et al.l 
(|2004l ) . Having specified the total dark matter dis- 
tribution, the caustic size, the magnification and 
the relation between halo mass and the image sep- 
aration need to be calculated numerically. We de- 
scribe the procedure we use, in the appendix A. 

O06 assumed the background source population to lie 
at a redshift of 2;^ = 2, and the source luminosity func- 
tion, ^( zg,I/) (X appropriate for the radio su rvey 
CLASS (|Mvers et aLl[2QQl iRusin fc Tegmarkll2QQl ) Q 
We are able to reproduce the expected ISDs from O06, 
given by Eq. [3 for all the three density distributions 
mentioned above. The expected ISDs corresponding to 
these three density distributions are shown in Fig. [9l 

3.4.2. Observed distribution 

The observed ISD is calculated by logarithmically bin- 
ning the image separations of 125 SARCS candidate^ 
with Ocut ^ 20 pixels (that is, > 3.7'') and an average 
ranking of 2 and above. The image separation for each 
lens candidate is taken as twice the Einstein radius or 
roughly the arc radius which is the distance between the 
candidate lensed image and the center of respective lens 
galaxy. Let — dO = 9\ and 0-\-dO = ^h, then the observed 

The luminosity density for such a steep faint end slope (q;lf < 
—2) diverges as L — )► and necessarily requires a cutoff below some 
value of I/min- 

Two of the candidates are excluded since their lensing con- 
figurations or the centers of their lens potential were ambiguous. 
These candidates have Ra = 0.0'' in Table E] 



ISD is given by 

^-^lens 



d9 



obs 



^ P(> ^h) - P{> 0i) 
N{> ^h) - N{> e,) 

N{> ecut)2d0 ■ 



(11) 



where the total number of observed lenses is A''(> 6'cut)- 
While comparing their theoretical predictions {dP/ dO) to 
the observed ISDs from the CLASS sample, O06 assumed 
an arbitrary normalization for their data points. Instead 
we note that, 



P\eU> 0) = 



iVlens(> 0) 
N\ens{> ^cut) 

Arie„s(> e) 



dP 



d0 



Pi> 0) 

P{> Ocut) ' 

dPi. 



Niens{> Ocut) 



d0 



P{> Ocut) ■ 



(12) 
(13) 



To facilitate a direct comparison of the observed ISD 
to the theoretical expectation from O06 (see Eq. [7]), we 
multiply Eq. [niby P(> 6>cut). The quantity P(> 6>cut) 
is obtained by integrating Eq. [3 from ^cut to 00. The 
SARCS data points are shown in magenta in Fig. [9l 
The vertical error-bars are calculated assuming Poisson 
number statistics and the horizontal bars show the bin 
width. 

The SARCS data points demonstrate that the average 
density profile of the halos, giving rise to the intermediate 
6 values (~ 3'^ — 12^'), is best represented by a combined 
profile for the main galaxy and the dark matter halo 
as opposed to a pure SIS or pure NFW profile. It is 
clear from the Fig. [9] that the SARCS sample follows a 
steeper ISD at the intermediate scales compared to halos 
with NFW profile and shallower than halos with only SIS 
profile. 

It is interesting to compare the SARCS sample with 
the SLACS and MACS samples which span the lower end 
and higher end of the ISD, respectively. These samples 
also have different selection functions compared to the 
SARCS sample. We apply the same procedure to calcu- 
late the observed ISD for SLACS and MACS data points 
shown in Fig. [9l Intriguingly, the SLACS data points lie 
in the regime where the SIS density profile just ceases 
to be dominant. Although the SLACS sample appears 
to be incomplete by a factor of 2 to 3 in the lowest 
bin, the ISD of SLACS could be seen as an extrapola- 
tion of the ISD from the SARCS sample. As expected, 
the MACS sample is nearly consistent with either the 
NFW or total profile. Assuming that incompleteness is 
the only major factor in the ISD of SLACS, the expected 
ISD corresponding to the total profile best matches the 
SLACS, SARCS and MACS samples combined. 

We note that we have not accounted for any effects due 
to purity or incompleteness of the SARCS sample in this 
paper. We suspect that the completeness of the sample 
as a function of the image separation is not severely af- 
fected due to the selection function. We are currently 
investigating this issue and the results will be presented 
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Fig. 10. — Image separation distribution. Theoretical curves for the total profile, shown with solid black line, are the same as in Fig. [9] 
Left: Adopting the c-M relation of MOT (dashed) steepens the total profile. The dotted curve, with M-L relation of Y08 and c-M relation 
of MOT, steepens further but negligibly. The total profile without the AC (dashed-dotted) is also shown. Right: The effect of varying q;lf 
on the ISD is shown for the total profile. The total profile without any bias (long dashed) is independent of q;lf- Changes in various model 
parameters have degenerate effects on the expected ISD. The current uncertainties on the data do not have further constraining power on 
the tested model parameters. All the theoretical curves are multiplied by the P(> 0^ ^n) since the SARCS sample consists of lenses with 

> 3.r^ 



in a forthcoming paper. 

3.4.3. Tests with varying models 

Here, we test the effects of varying the different com- 
ponents of the O06 model and compare against the ob- 
served ISD. Since different models have different values 
of P(> Ocut)^ we predict dP\Qns/dO from each model and 
compare it with the observed ISD in Fig. [TOl We test 
only for image separations spanning the observed val- 
ues. In the Fig. [9] and both the panels of Fig. [TOl the 
solid line represents the same ISD corresponding to the 
total two-component profile which accounts for the con- 
traction of the dark matter. First, we test the influence 
of excluding the AC while computing the total profile. 
This has the effect of making the expected ISD shallower 
as shown by the dashed-dotted curve in the left panel of 
Fig. [TOl Prima facie, the AC model fits the data bet- 
ter than the model without AC. However, as we show 
below, there are other degeneracies in the model which 
prevent us from ruling out the "no AC" model at high 
significance. 

Next, we test the effect of using different c-M rela- 
tions o n the ISD. For ex ample, we use the c-M rela- 
tion bv iMaccio et al.l (I2QQ7.. hence forth, MOT) instead of 
that given bv lBullock et al.l (|2QQl[ ). The c- M relation of 
MQ7 i s roughly 15-20% lower than that of Bullock et all 
(I2QQ1I ). The combined profile using the c-M relation of 
MOT is shown by the dashed curve in the left panel of 
Fig. [TOl Within the current statistical limits on the data, 
both the c-M relations appear plausible, although the 
data appears to slightly favor the c-M relation of MOT. 
Since the c-M relations differ significantly at small and 
large image separations, we will need additional samples 
of galaxy or cluster-scale lenses to test between the dif- 



ferent c-M relations. 

We also test the effect of using a more r ecent deter- 
minat ion of the M-L relation obtained by lYang etldl 
(|2008l . henceforth, Y08) from a sample of SDSS groups 
along with the c-M relation of MOT for the combined to- 
tal profile. The M-L relation of Y08 differ by -0.2 dex 
from that of I Vale fc OstrikeB ([2004) at the intermediate 
mass regime which is the regime of interest. This appears 
to cause a very negligible change in the predicted ISD. 

We try to quantify the effect of varying the slope of 
the source luminosity function at the faint-end. We show 
the effect on the combined profile of O06 and vary the 
power law index, of the source luminosity function, 
^(zs,I/) oc This influences the lens cross-section 

via the magnification bias. The solid curve in Fig. [TOl 
shows the expected ISD using the fiducial value of aLF = 
—2.1, while the dotted and short dashed curves show the 
ISD in the right panel, corresponding to aLF equal to -l.T 
and -1.2, respectively. It is evident from the figure that 
the observed ISD can be used to constrain the slope of 
the luminosity function, if the statistical error bars could 
be reduced. 

We note that the magnification bias factor in the bi- 
ased lens cross-section is calculated assuming that the 
background sources are point sources such as quasars. 
However, the background sources corresponding to the 
lensed arcs are mostly galaxies with extended surface 
brightness and their magnification bias could be negli- 
gible (e.g., Narayan & Wall ingtonl IT993) . Therefore, we 
calculate the ISD assuming no bias, that is, by substi- 
tuting B{zs) = 1 in Eq. [5l The long dashed curve in the 
right panel of Fig. [TOl shows the ISD without the bias. 

The data is consistent with all of the above tested mod- 
els within the uncertainties. The various scaling relations 
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from our model, that are tested here, have degenerate 
effects on the expected ISD. For example, varying the 
c-M relation has the same effect as changing the slope of 
the luminosity function or excluding AC. However, con- 
straints from independent observations such as dynamics 
and strong lensing could be used to determine the c-M 
relation. This will allow us to better constrain the slope 
of the luminosity function or make more robust state- 
ments about AC. 

For all the theoretical calculations, we have assumed 
z = 2 for the source redshift. We checked the effect of 
adopting different source redshift s, Zs = 1.5 and = 3.0 
on the predicted ISD. We found that using a higher 
(lower) source redshift causes the predicted dP/dO to 
be larger (smaller) by roughly a constant factor at all 
values of 0. However, the predicted dPiens/d^ is not af- 
fected because the corresponding increase (decrease) in 
^(> ^cut) almost perfectly cancels out the increase (de- 
crease) in dP/dO. This implies that the expected ISD 
would not be drastically different had we accounted for 
the distribution of source redshifts instead of assuming a 
single value for the source redshift. 

4. SUMMARY 

We have presented the SARCS sample from the com- 
pleted CFHTLS-WIDE and CFHTLS-DEEP covering 
a combined unmasked area of ~150 deg^ in the sky. The 
lens sample is compiled through a semi-automatic tech- 
nique consisting of using ARCFINDER algorithm, followed 
by visual inspection and ranking of the candidate s. We 
briefl y described the working of the arcfinder (jAlardI 
l2QQ6f ) and the modifications implemented in the newer 
version of the algorithm. Although the arcfinder V2.0 
is faster, there is still scope for improving the algorithm 
in terms of increasing the purity without compromising 
the completeness of the arc detections. 

We have compiled a total of 127 candidates in the 
SARCS sample, out of which 38 are found serendipi- 
tously. From the complete sample, 54 are promising lens 
systems (ranking of 3 and above). A total of 31 systems 
are almost certain or confirmed systems, out of which 
27 systems have been followed- up via techniques such as 
spectroscopy, high-resolution imaging and/or lens mass 
modeling. We found 2 radial arc candidates in the 
SARCS sample and both of them are located in the sys- 
tem SAIOO. The second radial arc is clearly identified 
from the high resolution HST imaging only. Our sam- 
ple may have more radial arcs which could be discovered 
with high resolution imaging. With statistics of radial 
and tangential arcs from a homogeneous and a larger 
sample of lenses, interesting constraints could be placed 
on the slope of inner density profiles of the dark matter 
distribution. 

We have discovered a total of 12 giant arcs (l/w > 8) 
in our sample. Statistics with giant arcs is considered 
to be a good probe of cluster properties or cosmological 
parameters but not an easy one. We have presented the 
redshift distributions of the lens galaxies with giant arcs 
and all arcs in the SARCS sample using the photometric 
redshifts and found to have mean values at z ~ 0.6. This 
is somew hat higher than the expec ted peak at redshift of 
0.3 — 0.4 (jBartelmann et aT]ll998f ) but consistent within 
the uncertainties. Note that the predictions need to be 
revised with improved simulations and more realistic as- 



sumptions. We also calculated the angular distribution 
of giant arcs which are sensitive to the ellipticity of the 
halo. We found an anisotropy in the orientation of arcs 
in our sample consistent with that seen in simulated lens 
clusters (e.g. see D04). In addition, the angular distribu- 
tion of all arcs from the SARCS sample exhibits similar 
anisotropy. This anisotropy appears to hold at all arc 
radii within the current uncertainties albeit needs to be 
verified with a sample of confirmed lenses. Thus, we have 
extended this result to group scale halos, observationally. 
It would be interesting to check if lensed arcs in simu- 
lated group scale halos also show a similar anisotropy in 
their azimuthal distribution and what constraints could 
be placed on the baryonic physics important in the inner 
few Kpcs which is probed by these lensed arcs. 

We followed the formalism of O06 to calculate the ex- 
pected ISD in order to compare it with the observed ISD 
from SARCS sample. We first reproduced the results of 
O06 and then introduced variations in different scaling 
relations used in their model. The SARCS sample al- 
lowed, for the first time, to probe the intermediate mass 
regime corresponding to group-scale halos via the ISD. 
We showed that the density profile of the halos are well- 
reproduced by a combined profile (NFW and Hernquist) 
at the group-scales, which is consistent with the predic- 
tions. Given the uncertainties in the data and the de- 
generacies in the model, both the models that account 
for or exclude the role of AC are consistent with the ob- 
served ISD. With the availability of larger statistics of 
confirmed lenses and understanding of the sample selec- 
tion function, the distinction between models with and 
without AC would be possible. 

Next, we varied the c-M relation, the halo mass- 
luminosity relation, the slope of a power-law source lumi- 
nosity function and the source redshift. We found that 
given the current uncertain ties in the obs erved ISD, the 
c-M re lations of both . Bullock" et al.l ()2001f ) and the more 
recent, [M accio et al.l (^20071 ). are plausible. The expected 
ISD does not vary significantly and fits the data well, if 
the halo mass-luminosity relations have an uncertainty 
of ~ 0.2 dex. Following O06, we adopted a power-law in- 
dex of aLF = —2.1 for the background source luminosity 
function to account for the magnification bias. We fur- 
ther tested the effect of varying the aLF on the expected 
ISD and found to be consistent with the data. However, 
within the uncertainties, the data is also consistent, if 
no magnification bias is assumed. We did not test the 
effects of any evolution of the luminosity function or any 
other functional form such as a broken power law. 

We found that varying the model parameters have de- 
generate effects on the ISD, for instance, changing the 
c-M relation and changing the slope of the luminosity 
function, aLF- Therefore, using priors from indepen- 
dent methods on one of these relations could help in 
constraining the others via the observed ISD. Since the 
background sources are assumed to be at a redshift of 2, 
we tested the effect of varying the source redshift. The 
expected ISD (dPiens/d6>) is not affected by choosing dif- 
ferent redshifts (zg = 1, 2, 3) between the range we tested. 

As described above, we have used arcs statistics to 
probe the average density profiles of group-scale lenses 
and we have shown the possibility to use arcs statistics in 
constraining some scaling relations. However, the models 
assumed in our work are simplistic and will need refine- 
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ment as the lens samples become larger with upcoming 
surveys. On the observational side, understanding the 
selection effects will also become crucial, if the model 
parameters need to be constrained with high accuracy. 
We hope to address some of these important issues in 
future studies. 
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APPENDIX 

A. NUMERICAL DETERMINATION OF LENSING QUANTITIES 

We describe how we numerically determined the relation between the halo mass and the image separation, the size 
of radial caustic and the position dependent magnification for general spherically symmetric mass distributions. These 
are used in the calculation of the expected ISD in the two of the models we tested, namely, the NFW and the total 
(NFW and Hernquist) profiles. 

Let the position of the lens in the plane of the sky be the origin of the coordinate system. The lensing equation 
relates the true angular position, /3, of the background source to the observed angular position of its image, in the 
plane of the sky, such that 

(3 = d-a{d), (Al) 

where a is the scaled defiection angle. In case of a spherically symmetric density distribution for the lens, the deflection 
angle a lies along the same direction as 6>, and the lens equation can be written as 

/3 = ^-^^, (A2) 

where the quantities in normal face are the magnitudes of the corresponding vectors in bold face. The relation between 
the deflection angle a and the magnitude of the image position is given by 

4GM(< Die) 

where Di is the angular diameter distance to the lens and M(< Di9) denotes the projected lens mass within a physical 
radius of size £^ = DiO^ which can be obtained using the following equation 

M(< = 27r j^^ R'dR' 2 p + R'^^ dz , (A4) 



where i^max = \^R^^^—^ and i?vir is the virial radius. 

In Fig. [TTJ we show the generic shape of how /3 and are related to each other for a spherically symmetric mass 
distribution which has a flnite density at the center (non-singular). For large values of 6 ^ /3 corresponding to the 
weak lensing regime. As the source approaches the lens in projection, i.e as /3 ^ 0, the source position /3 corresponds 
to multiple values of referred to as the strong lensing regime. The flrst instance of this is when the source is at the 
radial caustic (labeled as RC in Fig. [TT]) and here, the numerical derivative of P with respect to is zero. When /3 ~ 0, 
corresponds to the Einstein radius (Oe) which is the location of the tangential critical curve (see TCC in Fig. [TT]). 
The image separation is approximately equal to 20 e- 
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Fig. 11. — Source position as a function of image position for spherically symmetric lens density distribution. The center of the lens 
potential is at the origin. RC is the location of the radial caustic and TCC is the location of the tangential critical curve which corresponds 
to the Einstein radius. 



The magnification of the lensed image at d is given by the determinant of 



df3 



de 



In component notation, the required derivative is given by 

df3i , a\ d (a ^\ dO 



de 0^' 



30, 



The inverse of the determinant of the above equation simplifies to 



(A5) 

(A6) 
(A7) 
(A8) 

(A9) 



To summarize, given a spherically symmetric density distribution p(r), we first obtain the value of /3 for different 
values of 6 using equations IA2ilA4[ We use cubic splines to approximate the function /3{6). The value of /3 which 
corresponds to the local extrema of this function gives the angular size of the radial caustic. The image separation 
is given by two times the value of when /3 = 0. The mass within this angular image separation can be calculated 
using Equation IA4[ Finally, the magnification as a function of can be obtained using Equation IA9[ The required 
derivative is calculated by using the cubic spline fit. The magnification bias can thus be calculated using Equation [6l 

B. MEASUREMENT OF THE PEAK OF THE REDSHIFT DISTRIBUTION 

Consider a lensing galaxy with true redshift Zf for which the photometric redshift estimate is ± . Let us assume 
that the true redshifts of the lens population follow a Gaussian distribution, G(z^,cr|), with the mean and variance 
given by Zt and <j^, respectively. The probability that the estimated redshift of the lens galaxy is equal to Zi^ is then 
given by 



P{z,\G) = I P{z,\zt)P{zt\G)dzt, 

\Zi - Ztf 



1 

2^ 



-exp 



2a? 



{zt - Zt ) 
2a? 



2 \ 1 



dzt . 



Rewriting the factor inside the exponential as a quadratic in Zt and completing the square yields 

i2 



1 



271(7 i(7t 



exp 



-{Zj - ZtY 



exp 



1 



dzt . 



where z = (a| + a?) ^{zicr^ ^ + ZtCfi Carrying out the integral gives 



P{zi\G)-. 



zexp 



-{Zj - Ztf 

2K2 + <t2) 



(Bl) 

(B2) 
(B3) 

(B4) 
(B5) 
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Since the determinations of the photometric redshifts of lens galaxies are independent of each other, the combined 
likelihood for the data, given our model, can be written as jSf = Ylf=i P{zi\G). The posterior distribution for our 
model parameters given the data is then given by Bayes' theorem 

P{zua^\zi,a,)^^P{zuC7^) (B6) 

We assume uninformative priors on the parameters Zt and and sample from the above posterior distribution using 
a Monte Carlo Markov Chain. We quote the mean of the redshift distribution and the 68 percent confidence interval 
on it using the samples from the chain. 



